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FOREWORD 


Thls>  report  was  prepared  for  the  U.S.  Air  Force  Flight  Dynamics  Laboratory, 
Wri^t-Patterson  Air  Force  Base,  Ohio,  under  Contract  Number  F3361S-75-C-30S2. 

The  work  reported  herein  was  performed  at  the  Orlando  Division  of  Martin 
Marietta  Aerospace  as  a part  of  Project  8219,.  "Stability  and  Control  for  Aero- 
space Vehicles",  Work  Unit  82190117,  "Aerodynamic  Stability  Technology  for 
Maneuverable  Missiles". 

This  work  was  performed  during  the  period  February  1975  to  Deceiid>er  1976. 
The  principal  Investigators  were  J.  E.  Fldler  and  G.  F.  Aiello.  The  technical 
monitor  for  AFFDL  was  Kr.  WUllam  H.  Lane. 
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SUMMARY 


A computerized  engineering  model  is  presented  for  estimating  the 
effects  of  asymmetric  lee-side  vortices  on  slendet  missile  configurations. 

The  procedure  was  developed  using  both  empirically  determined  quantities 
and  t'  eoretlcal  techniques.  Empirical  inputs  define  both  vortex  locations  and 
street  vortex  strengths;  whereas,  potential  flow  considerations  guide 
in  the  definition  of  initial  vortex  strengths  and  Induced  forces  and 
moments.  The  procedure  Is  applicable  to  bodies  with  and  «dthout  tails. 
Calculable  effects  are;  induced  side  forces,  yawing  moments,  tail  forces, 
and  rolling  moments.  The  procedure  was  applied  to  a number  of  different 
ccr.b'i nations  of  geometries  and  flow  conditions  and  the  results  compared 
against  experimental  data.  These  comparisons,  while  not  exact,  have 
shown  the  procedure  to  be  suitably  accurate  for  preliminary  design  purposes. 
Using  this  procedure,  a user  can  estimate  the  magnitude  but  not  necessarily 
the  direction  of  vortex  induced  forces  and  moments  and  the  angle  of  attack 
at  which  they  first  appear.  Uncertalnlty  in  direction  is  attributed  to 
the  randomness  associated  with  formation  and  subsequent  shedding  of  the 
Initial  pair  of  vortices.  Nose  geometry  irregularities  greatly  Influence 
the  side  of  the  body  from  which  the  Initial  vortex  separates. 

A users  manual  for  this  computerized  procedure  is  also  provided.  The 
manual  Includes  user  Instructions,  a program  listing,  sample  Inputs, and 
sample  outputs. 
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1.0  INTRODUCTION 

In  recent  years,  increased  maneuverability  requirements  for  air-to-air 
■Issiles  have  dictated  corresponding  increases  in  angles  of  attack,  partlcur  , 
Icrly  for  those  vehicles  which  perform  slewing  maneuvers.  Maximum  angles 
for  such  missiles  can  now  reach  180  degrees.  As  a result  of  these  develop- 
ments, the  flow  fields  with  which  the  aerodynamlcist  has  to  deal  are  much 
■ore  complex,  and  traditional  methods  for  predicting  aerodynamic  character- 
istics are  Inadequate  to  deal  with  all  the  problems  involved.  One  example 
of  method  deficiency  is  found  in  the  angle  of  attack  range  25-50  degrees. 

Bere,  steady,  asymmetric  vortex  patterns  usually  develop  in  the  body  wake 
and  can  induce  large  side  forces  and  yawing  and  rolling  moments  which  are 
detrimental  to  missile  controllability.  Few  methods  are  currently  available 
for  predicting  these  forces  and  moments  or  for  analyzing  the  complex  flow 
fields  which  produce  them.  The  effects  are  particularly  acute  in  the  aub- 
sanlc-trsnsonic  speed  range.  At  Mach  numbers  greater  than  about  1.5,  however, 
they  tend  to  diminish  rapidly. 

The  flow  phenomena  involved  and  their  effects  have  received  considerable 
(1-9) 

attention  recently,  although  wake  vortex  effects  and  asynnetry  have 

been  reported  by  earlier  authors. Among  the  earliest  studies  were 

those  of  Allen  and  Perkins  and  Perkins  and  Jotgensen, in  which  some 

of  the  basic  flow  structure  in  the  wake  was  investigated  and  the  onset  of 

vortex  asymmetry  observed.  More  recent  work  by  Thomson  and  Morrison ^^^and 
(2  8) 

Thomson  ’ determined  details  of  the  wake  flow  field  and  the  associated 
vortex  characteristics  through  direct  flow  probing.  In  that  work,  attention 
was  drawn  to  the  strong  similarities  between  the  three-dimensional 
ssynietric  wake  and  its  two-dimensional  counterpart,  the  von  Karman  vortex  street 


1 


Measurements  were  made  of  vortex  strengths,  spacing  shedding  frequency.  The 
dats  of  Reference  1 have  often  formed  the  basis  for  follow-on  work  to  determine 
vortex  effects  on  slender  mlssilea^^*^^  or  for  the  computation  of  wake  flow 
characteristics. The  measurements  form  part  of  the  basis  for  the  present 
work,  with  suitable  modifications  for  various  flow  parameter  changes. 

Recent  wind  tunnel  investigations  performed  by  Martin  Marietta  showed 
marked  evidence  of  asymmetric  vortex  effects.  The  magnitudes  of  typical 
Induced  quantities  are  shown  in  Figures  1 and  2.  In  Figure  1,  pitch  and 
yaw  plane  moment  data  up  to  60  degrees  angle  of  attack  are  shown  for  an 
isolated  body  composed  of  an  ogive- cylinder  of  10:1  total  slenderness  ratio. 
Beginning  around  25  degrees  angle  of  attack,  considerable  yawing  moments 
were  Induced  due  to  asymmetries.  The  addition  of  tails  to  the  body  (Figure 
2)  results  in  the  generation  of  yawing  moments  which  are  almost  as  large  as 
the  pitching  moments.  In  addition  r.o  the  problem  of  magnitudes  Is  that  of 
unpredicatable  sign.  Random  changes  in  direction  of  forces  and 

moments  have  been  observed,  sometimes  related  to  changing  flow  conditions, 
but  often  Interpreted  to  be  caused  by  small  manufacturing  imperfections  near 
the  missile  nose  ^^*^’^\ln  fact,  significant  changes  in  Induced  force  and 
moment  magnitudes  and  signs  can  be  produced  by  rotation  of  all  or  parts  of 
the  body.^^’^*^^  This  has  further,  serious  Implications  for  missile  con- 
trollability. 

It  is  clear  that  techniques  are  required  for  calculating  the  forces 

and  moments  Induced  by  asymmetric  vortex  wakes.  Some  work  has  already  been 
(8  9) 

done  in  this  area  * ' for  slender  bodies  at  subcrltical  crossflow  Reynolds 
number.  The  proceduie  presented  In  this  document  deals  with  that  case  also, 
but  goes  further  in  considering  supercritical  crossflow  Reynolds  and  Mach 
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number  effects  for  bodies  with  and  without  tails.  This  technique  is 
aenl-emplrical , drawing  upon  the  experimental  evidence  referred  to  above, 
but  modifying  it  for  different  flow  conditions  and  supplementing  it  with 
analytical  results  and  techniques.  To  make  this  technique  more  readily 
useable,  it  has  been  programmed  for  digital  computation. 

The  layout  of  this  document  la  as  follows:  first,  a general  descript- 
ion of  the  flowfield  is  presentsd,  followed  by  descriptions  of  the  techniques 
used  to  model  the  asymmetric  vortex  wake  and  calculate  the  induced  forces 
and  moments.  Following  these,  there  is  a section  presenting  corjari- 
sons  between  predicted  results  and  experimentally  measured  data.  Finally, 
there  la  the  information  necessary  to  operate  the  program,  l.e.,user 
instructions,  a program  flw  chart  and  listing,  sample  inputs/outputs  and 
an  indication  of  program  limitations. 
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2.0  FLOWFIELD  DESCRIPTION 

When  a alender  missile  body  Is  placed  at  angle  of  attack  In  a uniform 
flow,  the  boundary  layer  generally  separates  on  either  side  of  the  body  and 
forma  a lee  side  wake.  Separation  usually  begins  near  the  rear  when  the 
missile  reaches  about  6 degrees  angle  of  attack.  The  wake  takes  the  form 
of  a pair  of  symmetrically-disposed,  counter-rotating  vortices  fed  by  vor- 
ticlty  shed  from  the  separating  boundary  layer.  As  angle  of  attack  increases 
the  axial  extents,  sizes  and  strengths  of  the  vortex  increase  also. 

In  general,  vortex  size  and  strength  also  Increase  towards  the  rear 
of  the  body.  Several  authors  have  formulated  descriptions  of  vortex  develop- 
ment along  slender  bodies  in  terms  of  two-dimensional,  impulsively-started 
flows  around  cylinders.  These  formulations  relate  flow  development  with 
time,  measured  either  from  the  beginning  of  impulsive  two-dimensional 
motion  or  from  the  instant  a fluid  particle  makes  contact  with  a three- 
dimensional  body.  In  the  latter  case,  time  is  defined  by  distance  travelled 
along  the  body  and  the  axial  component  of  freestream  velocity.  For  the  two- 
diiwnslonal  case,  th..  motion  of  the  vortex  cores  as  time  passes  (i.e., 
vortex  size  and  strength  increase)  theoretically  follows  a path  known  as 
the  FSppl  line.^^^^  Use 'has  been  made  of  this  result  in  the  present  work, 
as  will  be  described  later. 

When  the  body  angle  of  attack  reaches  about  25  degrees,  the  sytnnetrlc 
nature  of  the  wake  disappears.  The  two  vortices  are  Joined  by  a third, 
beginning  again  at  the  body  rear,  and  the  wake  becomes  asymmetric.  As 
angle  is  Increased  further,  more  vortices  Join  the  flow,  until  the  wake  con- 
tains several  which  have  been  shed  from  the  body.  An  idealized  model  of 
the  flow  field  is  shown  in  Figure  3.  A section  taken  through  the  wake  shows 
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VORTEX  CORES 


It  to  rosombla  the  von  Karmen  vortex  street,  trail  kk.(''yr  in  the  literature 
on  two-dimensional  flows.  Not  all  nf  the  vort  It'os  t)ri<  fully  alied,  liuwevcr; 
two  usually  remain  close  to  the  body,  receiving  vorticity  from  Che  sheddit^g 
boundary  layers.  If  angle  of  attack  continues  to  increase  these,  too,  will 
be  fully  shed  into  the  wake  when  their  strength  reaches  some  critical  value. 
Their  places  will  be  taken  by  yet  other  growing  vortices. 

While  a growing  vortex  la  receiving  vorticity  from  Che  boundary  layer, 
l.e.,  before  it  reaches  full  strength,  it  tends  Co  stay  close  to  the  body. 
Thomson  and  Morrison showed  schlieren  photographs  which  Impllad  that  the 
vo^rtex  core  is  somewhat  curved  while  it  is  forming.  Not  until  full  strength 
has  been  reached  does  Che  core  become  straight  and  parallel  to  Che  rest  of 
Che  vocticea  in  the  street.  This  core  behavior  is  indicated  in  Figure  3. 

It  has  been  determined  that  a marked  aimilarlty  exists  between  the 
auuiifeaCations  and  effects  of  two-  and  three-dlDenslonal  asymmecrlc  vortex 
wakes. In  fact,  Che  von  Karman  street  stability  criterion  that  the  ratio 
of  lateral  street  dimension,  h,  to  the  distance  between  vortices  of  like 
sign,  given  by  h ■ 0,281t,  has  been  shown  to  apply  to  both  cases 

(Figure  3).  However,  in  the  wake  of  Figure  3 no  lateral  motion  of  the  vor- 
tex cores  relative  to  the  body  t.ikes  place.  The  two-dimensional  phenom-.; ton 
of  increasing  distance  between  vortex  core  and  body  as  time  increases  is 
analoglsed  in  Che  same  way  as  for  the  symmetrical  vortices.  That  is  to  ssy, 

Che  wake  is  steady  and  the  motion  of  a fluid  particle  along  a (stationary) 
vortex  core  may  again  be  described  in  equivalent  time  by  its  axial  velocity 
end  distance  traveled.  Thrrugh  use  of  this  analogy,  Thomson^^^  and  Thomson 

and  Morrison  were  able  to  deduce  Che  strengths  of  asymmetric  wake  vortices 

(8) 

as  well  as  their  effects  upon  body  crossflow  drag.  Thomson  ^ ' recently  ex- 
tended this  work  Co  deal  with  induced  side  forces  and  yawing  moments  on  bodies 
at  subcrltical  crossflow  Reynolds  numbers. 


At  still  higher  angles  of  attark,  say  greater  than  50  degrees,  the  wake 
begins  to  display  unsteadiness.  The  vortex  cores  show  definite  lateral  dis- 
placeaents  relative  to  the  body  and  tlio  induced  forces  and  ooments  become 
tiM'xiependent . The  procedure  presented  in  this  document  does  not  consider 
this  case.  Attention  is  directed  only  to  the  phenomena  in  the  angle  of 
attack  range  from  25  to  SO  degrees  and  their  effects  upon  missile  aero- 
dynamic characteristics. 
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3.0  CONSTRUCTION  OF  FLOW  FIELD  MODEL 

Development  of  a procedure  to  predict  Che  effects  of  an  asymmetric 
vortex  wake  on  configuration  aerodynamics  requires  the  development  of  a 
realistic  flow  field  model.  A flow  field  model,  reflecting  the  state-of-the 
art,  can  be  developed  using  seml-emplrlcal  inputs  based  on  Che  data  of 

References  1,  2,  8 and  16  and  theoretical  results 

Vortex  strengths  and  locations  have  been  made  compatible 

with  the  findings  of  Reference  1 and  suitably  scaled  to  broaden  their  ranges 
of  applicability.  In  addition,  theoretical  results  on  the  locations  of  wake 
vortices  and  their  Images  have  been  introduced  and  Che  contribution  of  nose 
potential  lift  to  vortex  strength  has  been  considered.  Both  shed  and  grow- 
ing vortices  are  treated,  making  use  in  part,  of  a vortlclty-conservatlon 
concept. Each  of  the  above  components  of  the  model  is  described  in  de- 
tail beginning  with  vortex  strength. 

Vortex  Strength 

(1  8) 

Detailed  flow  surveys  have  shown'**  that  not  all  the  wake  vortices 
are  of  Che  same  strength.  Generally,  the  first  vortex  from  the  nose  origin- 
ates near  the  nose/body  Junction  and  has  the  smallest  strength  in  the  wake, 
r^,  say.  The  second  vortex  separates  soon  after  Che  first  and  has  a somewhat 
higher  strength,  r2.  From  the  third  vortex  onwards,  all  have  approximately 
the  same  strength,  r^(>r2),  a.id  their  spacing  and  strength  are  analogous  to 
those  of  the  vortices  in  a von  Karman  street.  While  the  first  and  second 
vortex  strengths  will  contain  contributions  from  the  (potential)  nose  lift, 
the  "street"  vortices  are  wholly  fed  with  vortlclty  from  the  separating  body 
boundary  layer.  Reference  8 presents  detailed  information  on  dimensionless 
vortex  strength,  for  various  angles  of  attack  at  subcritlcal  crossflow 
Reynolds  number.  For  the  first  two  vortices,  strength  is  calculated  using 
concepts  from  potential  flow  theory. 
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Th«  first  and  second  vortex  strengths  must  contain  a contribution 
froB  the  potential  flow  lift  of  the  nose.  The  Atrength  of  the  first  vortex, 
r^,  uy  be  eatlmated  If  It  Is  assumed  chat  the  nose  Is  replaced  by  a horizon- 
tal lifting  line  of  constant  strength,  located  at  the  nose/body  Junction 
(Figure  4).  All  of  the  nose  lift  la  assumed  generated  by  this  line  and  the 
associated  trailing  vortices  will  have  lea  strength  provided  th^’y  receive  no 
further  vorticlty,  from  the  shedding  boundary  layer  for  example.  The  first 
vortex  shed,  near  the  nose  body  junction,  will  probably  contain  only  potential- 
flow-generated  circulation.  To  calculate  vortex  strength  is  straight 
forward.  It  can  be  shown that  the  coefficients  of  normal  and  axial  force 

acting  on  the  nose  are  predicted  by  slender  body  theory  to  be; 

2 

Cjj  ■ 2 sin  a and  “ sin  a 

o 

where  the  reference  area  la  that  of  the  body  cross  section.  Converting  these 

quantities  to  lift  coefficient  yields 

C,  ■ C„  cos  a ~ C.  sin  a 
L ft  A . 

If  C.  is  assumed  negligible  (as  it  will  be,  compared  to  the  axial  forces 
*o 

generated  at  the  high  angles  of  attack  here)  this  expression  becomes  after 
some  manipulation 

L " ® sln^  a) 

- P r^dv 

The  latter  expression  Is  the  well-known  Kutta-Joukowskl  theorem.  Finally, 
the  equation  may  be  rearranged  to  yield: 


Vd  sin  a 


« 2 
■ ^ (2  cos  a + sin  a) 


(1) 


This  dimensionless  vortex  strength  may  be  compared  with  measurements 
from  flow  field  surveys. Using  an  angle  of  30  degrees,  equation  (1) 
yields  a dimenslonlesa  strength  of  0.78.  The  corresponding  measured  strengths 
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from  Reference  1 show  a range  of  values  from  0.3  to  1.1  with  a mean  value 
of  0.7.  This  compares  quite  well  with  the  predicted  value.  Where  it  is 
found  that  the  first  vortex  is  shed  forward  of  the  nose/body  Junction,  its 
strength  is  simply  determined  by  using  the  local  nose  diameter  in  equation 


There  should  be  little  or  no  effect  of  Reynolds  number  on  since  the 
strength  is  determined  by  potential  flow  considerations.  On  the  other  hand, 
is  wholly  produced  by  vlsous  flow  and  is  strongly  affected  by  crossflow 
Reynolds  number  (R^),  which  influences  the  characteristics  of  the  boundary 
layer  shed  to  form  the  street  vortices.  The  second  vortex,  being  shed 
downstream  from  the  nose/body  junction,  is  modeled  here  as  a ''mixed  vortex," 
i.e.,  a potential  flow  vortex  which  receives  additional  vorticlty  from  the 
shedding  boundary  layer  just  aft  of  the  nose  (Figure  4).  Formally,  is 
expressed  as: 

r^-Ti  + Kr^ 

K tras  found  empirically  to  be  about  0.22  from  the  wake  survey  data  of  Ref- 
erence 1. 

For  street  vortex  strength  at  subcritical  , the  data  shown  in  Figure 
(8 ) ^ 

5 were  used  ' In  order  to  scale  F for  supercritical  R , use  was  made  of 

® © 

(14)  c 

von  brman's  result  that  crossflow  drag  coefficient,  is  approximately 

proportional  to  the  street  vortex  strength,  F (the  expression  for  C.  con- 

9 OC 

2 

tains  terms  in  F^  and  F^  ; however,  the  latter  accounts  for  less  than  10  per- 
cent of  the  total;  hence,  is  approximately  proportional  to  F^).  Data  obtained 
from  Martin  Marietta  investigations  into  crossflow  drag^^^^  produced  the  in- 
formation of  Figure  6,  which  shows  that  for  low  crossflow  Mach  number,  crossflow 
Reynolds  number  has  a strong  Influence  on  crossflow  drag.  Above  the  critical 
crossflow  Reynolds  number  (about  10^),  shows  a significant  decrease  below 
the  subcritical  value.  Hence,  F^,  too,  will  be  significantly  reduced.  An 

Increase  in  crossflow  Mach  number,  M , is  required  to  Increase  bbth  C.  and 

c ac 


' ' ' 
i--  .1  . - ' * . 


r . Accordingly,  for  R > 10  the  sober  ideal  R values  of  T from 
* *c  *c 

Figure  5 are  scaled  in  the  same  ratio  as  from  Figure  6. 

In  order  to  illuminate  further  the  mechanisms  underlying  the  relation- 
ship between  crossflow  Reynolds  number  and  vortex  strength,  recent  work  by 
Fidler^^^  will  be  briefly  described.  It  nay  be  shown  that  the  strength  of 
a lee  side  asywetric  vortex  can  be  related  to  the  flux  of  crossflow  vorticity 
leaving  the  body  at  the  circumferential  separation  point  defined  by  0^.  Further, 
the  vorticity  flux  shed  over  axial  distance  g (at  fixed  angle  of  attack)  in 
unit  time  is: 


F,  • const  In  dx 

B / s 


The  vorticity  thus  shed  diffuses  through  the  vortex  to  produce  the  cir- 
culation r . Equating  the  flux  F to  that  flowing  along  a wake  vortex  of 

8 D 

core  velocity  yields  an  expression  for  vortex  strength: 


const 

V 


K 

/“■ 


Now  U^,'  the  circumferential  velocity  at  the  boundary  layer  edge  at  sepa- 


ration, is  a function  of  0 . In  the  two-dimensional  case 

a 


it  has  been 


found  experimentally  that  when  crossflow  Reynolds  number  is  subcrltlcal,  9^ 
lies  near  the  meridian  of  the  cylinder  and  the  associated  separation  velo- 
city is  larger  than  in  the  supercritical  Reynolds  number  case,  where  the 
separation  point  lies  far  over  bn  the  cylinder  lee  side.  Continuing  the 
anelogy  between  two  and  three  dimensions  indicates  that  crossflow  Reynolds 
number  affects  through  its  effects;  first  upon  9^,  which  in  turn  deter- 
mines which  defines  the  vorticity  flux  flowing  from  the  body  and  dif- 
fuses through  the  vortex  to  produce  P . Hence,  subcrltlcal  R produces 

c 

larger  P values  than  does  supercritical  R , provided  crossflow  Mach  number 

c 

is  low. 
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The  foregoing  discussions  deal  with  fully-developed  vo  ..i.ces  which 
have  left  the  Immediate  vicinity  of  the  body  and  Joined  the  asymmetric  wake 
pattern.  However,  consideration  must  be  given  to  the  street  vortices  while 
they  are  forming  close  to  the  body.  As  previously  discussed,  these  growing 
vortlclcs  are  fed  by  vorticlty  from  the  separating  body  boundary  layer. 
Assuming  rapid  diffusion  of  vorticlty  the  local  vortex  strength  r(x)  formed 
by  shedding  a boundary  layer  over  distance  x may  be  written  (see  previous 
equation) . 


r(x)  -const 


x 

/“ 


2 
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Furthermore,  it  has  been  experimentally  determined  ' that  for  those 
portions  of  the  body  where  street  vortices  are  being  formed  and  for  given 
flow  conditions,  the  circumferential  location  of  boundary  layer  separation, 

6 , is  approximately  constant  with  x.  Hence,  U , the  circumferential  ' 

8 8 

velocity  at  the  boundary  layer  edge  at  separation  is  also  constant.  Thus 
the  local  vortex  strength  is  seen  to  be  directly  proportional  to  x,  l.e., 
r(x)  “ x/g  where  g Is  the  distance  over  which  boundary  layer 
fluid  is  shed  to  form  a street  vortex. 

One  final  item  affecting  the  strengths  of  forming  vortices  has  been 

(8) 

considered.  It  has  been  shown  by  Thomson  that  as  a growing  vortex  nears 
the  base  of  the  body,  the  rate  at  which  Its  strength  Increases  Is  reduced 
due  to  base  proximity.  Thomson's  data  Indicate  that  the  growth  rate  Is 
only  ^0  percent  of  normal.  The  axial  extent  of  the  region  over  which  the  base 
Influence  Is  felt  Is  given  by 


2 S 1 tan  a 

The  model  contains  this  feature. 

At  any  body  axial  station,  then,  the  strengths  of  all  the  vortices 
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are  calculable  from  the  combined  theoretical/empirical  procedures  described 
above.  The  next  stage  In  flow  field  model  construction  Is  to  locate  the 
vortices  relative  to  the  body  so  that  their  effects  may  be  calculated. 
Vortex  Location  and  Spacing 

The  problem  of  locating  vortex  cores  In  space  relative  to  the  body 
must  be  handled  separately  for  vortices  growing  near  the  body  (i.e.,  those 
being  f^  with  vortlclty  from  the  separating  boundary  layer)  and  for  shed 
vortices  which  can  be  considered  part  of  the  wake  street.  In  the  former 
case,  use  Is  made  of  theoretical  results  from  ttio-dlmensional  potential 
flow  theory;  In  the  latter,  systematic  experimental  evidence,  suitably 
scaled  for  flow  parameter  changes.  Is  employed.  The  case  of  growing  vor- 
tices will  be  described  first. 

Growing  Vortices 

la  order  to  model  the  trajectory  followed  by  a growing  vortex,  two 
reference  points  are  required.  The  first  Is  the  clrcurn'erentlal  location 
on  the  body  at  which  the  boundary  layer  separates.  The  second  is  the 
point  In  space  at  which  the  vortex  has  reached  full  strength  and  can  be 
said  CO  have  Joined  the  street. 

The  first  point  is  located  by  the  empirical  relationship 

0^  ■ sin  ^ (3  tan  6/2  tan  o) 

The  angle  thus  defined  places  separation  in  the  general  region  typical 
of  laminar /turbulent  boundary  layers.  Since  the  means  of  determining 
vortex  strength  does  not  rely  upon  an  exact  knowledge  of  the  separation 
point  location,  the  above  approximation  is  sufficient  for  present  purposes. 
If  the  determination  of  vortex  strength  had  required  use  of  the  Kutta 
condition,  or  some  estimate  of  vortlclty  flux  leaving  the  body,  the  angle 


^ -a-**  • , - 


6^  would  have  been  required  with  some  accuracy.  In  Che  present  method 
however,  vortex  strength  la  defined  otherwise  and  a rough  estimate  of  the 
separation  angle  is  sufficient. 

The  second  point  for  anchoring  Che  growing  vortex  trajectories  was 
taken  as  the  Intersection  of  the  Foppl  and  von  Kerman  lines  along  which 
the  symmetric  and  asymmetric  vortex  cores  were  kno«m  to  move  respectively. 
Use  of  this  point  was  Justified  as  follows: 

As  discussed  earlier,  the  appearance  of  vortices  in  the  wake  shows, 
at  the  earliest  stage,  a symmetric  pair.  At  any  axial  station,  increasing 
angle  of  attack  produces  Increased  vortex  size  and  strength  as  well  as  an 
outward  movement  of  the  core  approximately  folloi^ing  the  Foppl  line.  It 
was  reasoned  Chat  at  Che  first  appearance  of  asymmetry,  one  of  the  vortices 
on  the  Foppl  line  «K>uld  change  its  trajectory  and  proceed  outwards  along 
a new  path  defined  by  the  von  Kaman  stability  relationship  described 
earlier.  The  second  vortex  would  perform  similarly  and  this  would  then 
set  up  the  spacing  of  vortices  in  the  street.  It  was  hypothesized  then, 
that  in  order  to  continue  the  spacing  pattern,  the  Intersection  of  the 
Foppl  and  von  Karman  lines  would  denote  the  point  at  which  all  vortices 
reached  full  strength  and  were  shed  into  the  street.  In  order  to  test 
this  hypothesis  the  schlleren  photographs  of  Thomson  and  Morrison 
were  examined  to  determine  the  distance  from  the  body  where  vortex  feed- 
ing from  the  boundary  layer  ceased,  i.c.,  the  point  at  which  the  growing 
vortex  cores  became  straight  and  Joined  the  street.  It  was  found  that  the 
points  thus  defined  covered  a band  of  values  from  2 to  2.S  diameters  above 
the  body  and  at  a lateral  distance  defined  by  the  von  Karman  stability 

criterion.  To  determine  the  theoretical  location  of  the  shedding  point, 

2 2 

the  following  procedure  was  used.  With  the  equation,  (14),  2 ry  ■ r - a , 
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Che  Foppl  line  was  drawn  relative  to  the  cylinder  of  Figure  7.  The  data 
of  Reference  1 were  used  in  Che  equation  li/l  “ S/x  Co  determine  t,  Che 
apaclng  between  street  vortices  of  like  sign  for  aubcricical  crossflow 
Reynolds  numbers.  Then  the  von  Karman  relation  h >0.2811  was  used  to 
superimpose  the  street  vortex  location  line  on  Figure  7.  The  Intersection 
of  Che  Foppl  and  von  Karman  lines  was  found  to  lie  within  the  experimentally' 
determined  range  given  above.  It  was  concluded  then,  that  this  intersection 
point  provided  a good  estimate  of  the  location  at  which  the  vortices  stopped 
growing  and  were  shed  to  form  part  of  the  street.  For  the  purposes  oi  this 
engineering  flow  modtl  the  growing  vortex  core  was  assumed  to  move  linearly 
between  the  two  anchor  points.  The  resulting  core  trajectory  model  is 
shown  (foreshortened)  in  Figure  8. 

Street  Vortices 

From  the  Foppl/von  Karman  line  Intersection  the  vortex  cores  stream 
back  into  Che  wake  as  straight  lines,  making  angle  ^ with  the  body  axis. 

This  angle  was  measured  by  Thomson  and  Morrison, related  to  the  rate  at 
which  a fluid  particle  flowing  along  the  core  increases  its  distance  from 
the  body,  and  thence  analogized  to  the  two-dimensional  von  Karman  street. 

( is  related  to  missile  angle  of  attack  through  the  parameter  x * tan 
'1) 

C/tan  a . If  the  point  on  the  body  from  which  a vortex  emanates  is 
known,  then  Its  core  location  In  space  may  be  partly  determined  using  x* 

For  the  purposes  of  this  work  x vas  assumed  unaffected  by  changes  in  cross- 
flow  Reynolds  number. 

Vortex  starting  points  on  the  body  were  estimated  in  Reference  1 as 
being  the  intersection  point  of  the  body  axis  and  the  extrapolated  street 
vortex  cores.  Since  these  positions  were  obtained  for  tailless  bodies,  the 
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FIGURE  7 COMPARISON  BETWEEN  PREOCTED  AND 

experimental  point  where  vortex 

IS  SHEO  INTO  STREET 
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FIGURE  8 CROWING  VORTEX's  TRAJECTORY 


present  wrk  eseunes  chat  In  the  presence  of  tails,  body  vortex  formation 
ceases  at  the  body/tall  leading  edge  Intersection.  Provided  the  Scrouhal 
Buaber,  S,  is  known,  Che  starting  points  of  wake  vortices  can  be  detsmined. 
Further,  using  the  relatlon^^^  d/1  • S/x,  the  lateral  spacing,  1,  between 
vortices  of  like  sign,  and  hence  0.2811,  the  von  |Carman-deduced  criterion 
for  lateral  spacing  can  be  determined.  The  vortex  cores  can  now  be  located 
in  space  relative  to  the  body  at  any  point  along  their  lengths.  At  this 
stage,  however,  the  Information  Is  sufficient  only  for  subcrltlcal  cross- 
flow  Reynolds  numbers. 

The  quantity  on  which  attention  must  be  concentrated  when  scaling 
vortex  locations  for  crossflow  Reynolds  number  Is  Che  Strouhal  number,  S. 

This  measure  of  Che  rate  at  which  vortices  of  like  sign  are  shed  from  Che 
body  has  been  shown  to  exhibit  strong  similarities  between  Che  two-  and  three 
dimensional  cases.  In  Reference  1,  S was  determined  experimentally  to  have 
values  near  0.2  for  a wide  range  of  crossflow  Mach  numbers.  This  compares 
well  with  Che  two-dimensional  value  at  subcrltlcal  crossflow  Reynolds  n^ber. 
It  Is  known  however,'  that  tor  cylinders,  S varies  with  crossflow  Reynolds 
number  and  increases  to  values  in  the  range  0.3  to  0.5  at  supercritical  Rec* 
This  has  a direct  effect  on  vortex  spacing,  both  longitudinal  and  lateral. 

For  bodies  alone  then,  S • 0.2  is  used  for  subcrltlcal  Re^*  been 

found  that  a value  of  S ■ 0.35  gives  Che  best  results  for  supercritical  Re,.* 
For  bodies  with  tails,  Che  race  of  vortex  formation  Is  Influenced  by  the 
presence  of  the  tails.  Experimental  evidence  Indicates  that  the  effect  is 
Co  reduce  the  rate  of  formation  to  chat  defined  by  S > 0.2  regardless  of 
whether  the  crossflow  Reynolds  number  Is  sub-  or  supercritical.  Accordingly, 
for  bodies  with  tails,  S <•  0.2  la  used  throughout. 
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Using  the  shove  enplrlcal  and  theoretical  inputs,  a model  can  be 
constructed  of  the  asymmetric  wske  produced  by  a slender  missile  configuration 
at  high  angles  of  attack.  This  model  Includes  the  number  and  locations  of 
vortices  in  the  wake  and  their  strengths,  suitably  scaled  for  Mach  number 
aa^  Reynolds  number  effects.  Both  gro%d.ng  and  shed  vortices  are  included 
in  the  model.  Having  constructed  a model  of  the  wake,  the  next  step  is  to 
consider  its  effect  on  the  missile  configuration. 
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4.0  FORCE  AND  MOMENT  COMPUTATIONAL  TECHNIQUES 

Having  astabllshed  a technique  to  model  the  wake  produced  by  a slender 
■Issile  configuration  at  high  angles  of  attack,  the  next  step  i>i  to  calculate 
the  forces  and  moments  Induced  on  the  configuration  by  the  wake.  This 
section  preoents  a computational  proc>  lure  to  do  this  for  bodies  alone  and 
then  for  bodies  with  tails.  In  either  case,  the  procedure  calls  for  the 
calculation  of  Incremental  effects  produced  on  various  configuration  segments 
and  then  integration  of  these  effects  to  determine  the  forces  and  moments 
induced  by  the  wake. 

This  computational  procedure  was  designed  with  the  intent  that  it 

not  be  overly  complicated..  With  this  in  mind,  the  following  basic  assumption 

«ias  made  concerning  the  location  within  the  body  of  the  image  vortices 

required  to  preserve  the  velocity  tangency  condition  on  the  surface  of  the 

body  at  each  axial  station.  These  image  vortices  are  located  in  the  body 

on  a plane  perpendicular  to  the  vortex  core;  however,  this  raises  a problem 

because  of  the  inclination  of  the  vortex  cores  relative  to  the  body  axis. 

A cross-section  normal  to  the  street  cores  shows  an  elliptical  body  section, 

inside  which  the  location  of  image  vortices  is  not  simply  accomplished.  In 

keeping  with  the  simple  nature  of  this  model  it  was  decided  that,  if  possible, 

(14) 

image  vortices  should  be  located  by  means  of  the  circle  theorem.  This 

was  accomplished  by  resolving  Che  vortex  circulation  vectors  normal  and 
parallel  to  Che  body  axis.  By  ignoring  the  former  as  having  no  relevance 
in  the  two-dimensional  section  model,  the  latter  components  plus  their 
Images  could  then  be  used  to  determine  forces  and  moments.  The  use  of  Images 
was  not  necessary  for  body  quantities,  but  was,  however,  mandatory  for 
tailed  regions  of  the  missiles.  These  points  are  discussed  below. 
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Body  Forces  and  Moments 

The  calculation  of  body  forces  and  moments  begins  with  the  determination 
of  the  Incremental  force  on  each  body  segment.  In  order  to  calculate  this 
force,  the  local  net  circulation  of  all  the  vortices  in  the  wake  must  be 
known.  Further,  these  vortices  must  have  their  circulation  vectors  resolved 
parallel  to  the  body  as  described  above.  At  any  axial  station  from  the  nose, 
there  are  usually  several  wake  vortices  of  strength  and  two  growing 
vortices,  one  on  the  left  side,  having  strength  rj^(x)  and  another  on  the 
right,  of  strength  Since  g is  the  axial  distance  over  which  a vortex 

grows  to  full  strength  then  x^^  - x^  -I*  g/2  assuming  To  resolve 

the  strengths  parallel  to  the  body  axis,  t,he  angle  a for  growing  and  C for 
street  vortices  were  used  respectively.  Side  force  on  unit  length  of  the 
body  was  calculated . using  the  Kutta-Joukowskl  expression. 

m R-r.  nrt-l  R-r 

6Y  - pVsln  a l(r-(x)-r,  (x)>  cos  <j  + { E (-r-^)r  - E } cos 

, 1-1  R “Rj  1-1  R *L^ 

Where  it  has  been  assumed,  for  illustration,  that  there  are  (m'fl)  fully 
developed  vortices  on  the  left  side  of  the  body  and  m on  the  right.  The 
term  (R-r^)/R  is  an  empirical  factor  accounting  for  the  attenuation  of 
fully  developed  vortex  effects  on  the  body  as  their  distances  from  it  increases 
r^  is  the  actual  distance  of  the  vortex  from  the  body  and  R is  an  arbitrary 
distance  at  which  vortex  effect  is  assumed  to  have  attenuated  to  aero.  The  ' 
aauilleat  value  of  r^  is  defined  by  the  Foppl/von  Karman  intersection  point. 

This  formulation  was  an  attempt  to  model  the  expected  vortex  attenuatlcn 
effect  empirically.  As  will  be  shown  later,  it  was  found  that  the  best 
representation  of  forces  and  moments  was  obtained  when  the  street  vortex 
strength  was  allowed  to  attenuate  to  zero  immediately  after  shedding. 
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Total  aide  force  Y and  yawing  moment  YM  are  ealculat eit  by  numerically 
integrating  along  the  body  length.  At  this  stage  In  the  compncerized 
version  of  the  procedure,  an  option  is  provided  allowing  the  user  to  cal- 
culate Che  effect  of  varying  nose  fineress  and  bluntness  ratios.  This 
option  uses  Che  data  of  References  3 and  19  to  scale  the  total  body  side 
force  and  yawing  moment  for  various  nose  fineness  and  bluntness  ratios. 

The  calculation  of  individual  side  forces  and  yawing  moments  outlined  above, 
employs  the  flowfleld  model  developed  urging  the  empirical  data  from  Reference 
to  determine  street  vortex  strength  and  vortex  separation  points.  These 
empirical  inputs  were  derived  from  wind  tunnel  data  for  a configuration 
having  a 3.798  fineness  ratio  conic  nose.  The  data  of  Reference  3,  obtained 
from  tests  of  configurations  with  noses  of  fineness  ratios  of  2,  3 and  4, 
show  chat  Che  maximum  absolute  value  of  side  force,  generally 

tended  to  increase  with  Increased  fineness  ratio  for  noses  with  little  or 
no  bluhCness  (see  Figure  9).  To  account  for  this  effect,  the  data  of  Figure  9 
were  used  to  produce  scale  factors  by  which  Che  basic  output,  based  on  a 
3.798  fineness  ratio  nose,  could  be  scaled. 

The  data  of  P.eference  3 also  indicated  that  increases  in  nose  bluntness 
generally  tended  to  red-ice  the  maximum  value  cf  side  force.  Figure  10 
shows  the  effects  of  blunting  a fineness  ratio  4.0  ncse.  Data  were  also 
available  for  fineness  ratio  2.0  and  3.0  nose  configurations  with  blunCnesses 
of  0,  5,  10,  20  and  50  percent.  These  data  were  used  to  produce  scale 
factors  by  which  the  values  of  side  force,  already  corrected  for  fineness 
ratio,  could  be  scaled. 

Schlieren  photographs  in  Reference  3 indicate  that  variations  in  fineness 
ratio  Impact  the  location  at  which  vortices  are  shed  and  that  variations  in 
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fineness  ratio  Impact  the  location  at  which  vortices  are  sh6d  and  that 
variations  in  nose  bluntness  delay  vortex  separation.  .This  is  to  say  that 
variations  in  nose  geometry  influence  the  structure  of  the  asymmetric  wake. 
Rather  than  attempting  («  modify  tlio  flowflold  model  presonlod  In  Suction 
3.0,  the  scaling  technique  outlined  above  was  employed  in  keeping  with 
the  simplified  nature  of  the  model. 

Caution  is  advised  when  employing  the  above  scaling  technique  for  nose 

fineness  an.d  bluntness  ratio  variations.  The  war.iing  is  necessary  because 
3 19 

the  data  * used  to  derive  the  scaling  factors  were  mean  values  and  data 
showing  considerable  scatter.  It  is  also  possible  that  an  Independent  test 
of  these  same  configurations  could  produce  considerably  different  data  due 
to  the  dynamic  nature  of  the  vortex  shedding  phenomenon. 

Tall  Forces  and  Moments 

Induced  side  forces,  yawing  moments  and  rolling  moments  due  to  the 

presence  of  tails  can  arise  even  when  no  asymmetry  exists  in  the  wake  vortex 

pattern.  This  is  because  of  varying  net  angles  of  attack  of  the  various  tails 

when  the  missile  is  rolled  at  arbitrary  angles.  Hence,  in  order  to  be  certain 

that  experimental  data  indicate  the  presence  of  vortex  asymmetry,  it  is 
necessary  to  be  selective  in  choosing  missile  roll  attitudes.  For  a cruciform 
missile,  roll  angles  of  0 "plus"  and  45  degrees  "cross"  are  the  only  attitudes 
where,  in  the  absence  of  wake  asymmetry,  sero  side  forces,  rolling  moments 
and  yawing  moments  occur  (assuming,  of  course,  no  tail  deflections).  The  pre- 
sent model  has  been  compared  against  experimental  data  from  cruciform  missiles 
and  hence  the  above  two  roll  angles  will  be  referred  to  exclusively. 

The  procedure  Cor  calculating  induced  tall  forces  will  be  described  for 
a vertical  lee-aide  tall  such  as  could  be  used  on  a cruciform  missile  in 
"plus"  attitude.  The  major  elements  of  the  treatment  are  most  expeditiously 
described  for  this  case,  although  the  complete  procedure  will  handle  tails  . 
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at  any  roll  angle  so  that  the  program  can  treat  ’'odles  at  arbitrary  roll  or 
whose  tails  number  other  than  four. 

The  first  step  in  the  process  is  to  determine,  for  each  vo*’tex,  the 
sldewash  velocity  induced  normal  to  the  tail  leading  edge  at  each  spanwlse 
station.  Since  the  effects  of  image  vortices  must  be  considered  and  it  is 
preferred  that  the  circle  theorem  be  used,  the  wake  vortex  circulation 
vector  is  resolved  parallel  to  the  body  axis  as  before.  In  this  case, 
however,  instead  of  simply  resolving  fg  through  some  angle  analogous  to  C 
as  was  done  for  the  body  ab6ve,  a new  stre.^gth  is  defined  which  will  pro- 
duce the  same  values  of  sldewash  velocity  as  did  the  original  vortex.  i 

A single  vortex  in  the  wake  plus  its  images  will  induce  some  distribution  i 

j 

of  sldewash  velocity  v along  the  tail  leading  edge.  The  average  value 
of  this  sldewash  is 

1 P ’ 
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The  total  Induced  forces  and  moments  are  obtained  by  summing  the 
effects  of  all  vortices  in  the  wake. 

Yawing  moment  induced  by  the  tall  is  calculated  by  assuming  that  the 
resultant  loading  is  located  at  the  mid-point  of  the  leading  edge. 

The  above  procedure  has  been  generalized'  to  handle  tails  at  any  roll 
angle.'  In  this  way,  the  contributions  of  all  calls  to  side  forces  and 
ya%rlng  and  rolling  moments  may  be  determined. 

Since  arbitrary  roll  angles  can  be  considered,  this  raises  potential 
problems  when  a vortex  core  intersects  the  leading  edge.  For  potential 
vortices  of  the  kind  used  here,  flow  velocities  become  extremely  high 
.near  the  core,  resulting  in  unrealistically  high  Induced  sidewash  angles. 
This  problem  was  circumvented  by  Introducing  into  each  vortex  a solid  core 
of  radius  0.25  body  radius. 

When  the  tall  Intersects  this  core,  calculations  are  discontinued.  If 
Che  core  passes  over  a Call  while  angle  of  attack  is  increasing,  the  forces 
and  moments  are  faired  across  the  gap  within  which  calculations  are  not 
performed.  This  procedure  has  proved  quite  satisfactory  in  practice. 

Instead  of  the  slender  body  theory  approach  to  tail  force,  strip 
theory  might  have  been  employed.  It  was  felt,  however,  that  the  former  was 
more  appropriate  for  Che  low  aspect  ratio  calls  for  which  the  program  would 
probably  be  used.  The  program  is  flexible  enough,  however,  that  strip 
theory  could  be  easily  introduced  if  it  were  considered  necessary. 
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5.0  COMPARISONS  BETWEEN  COMPUTED  AND  EXPERIMENTAL  RESULTS 


At  this  stage  then,  the  effects  of  the  wake  vortices  on  bodies  and 
tails  are  calculable.  Comparisons  will  now  be  presented  between  predicted 
and  experimental  forces  and  moments  to  illustrate  the  performance  of  this 
engineering  model.  Bodies  alone  and  with  cruciform  tails  will  be  consid- 
ered. The  effect  of  varying  the  number  of  vortices  considered  will  be  shown. 

Program  predictions  have  been  compared  against  various  experimental  data, 
generated  on  a selection  of  the  Martin  Marietta  Aerodynamic  Research  models. 

The  basic  model  used  for  high  angle  data  generation  was  a 10  caliber  tangent 
ogive/cylinder.  This  body  was  tested  alone  and  with  several  sets  of  cruci- 
form tails  affixed.  Each  tall  was  individually  Instrximented  so  that  program 

j 

predictions  of  vertical  tail  forces  and  moments  could  be  checked  for  "plus" 
attitudes.  For  tails  ocher  than  vertical,  the  program  would  predict  only 
the  incremental  force  generated  on  them  by  the  asymmetric  wake.  Freestream 
effects  were  not  taken  into  account,  and  hence  forces  for  non-vertical  tails 
will  not  be  discussed.  However,  since  the  differences  between  non-vertical 
tails  due  to  asymmetric  effects  were  presumed  to  be  valid  for  "plus"  and  "cross" 
attitudes,  the  program  was  used  to  predict  rolling  moments  for  these  attitudes. 

The  first  comparison  is  shown  in  Figure  11.  Here,  program  predictions 
of  normal  force  on  a vertical,  lee  side  tail  fixed  to  the  body  are  coBq>ared 
against  experimental  data.  Also  shown  is  Che  effect  of  considering  all  of 

i 

Che  vortices  in  the  wake  and  of  using  only  those  vortices  closest  to  the  1 

body,  l.e.,  the  growing  vortices.  It  will  be  seen  that  the  magnitude  of 
call  force  is  predicted  within  a few  percent  using  only  the  growing  vortices, 
while  use  of  all  vortices  produces  a significant  discrepancy.  The  angle  at 
which  asymmetry  begins  is  matched  only  fairly,  but  Che  angle  at  which  the 
appearance  of  a new  wake  vortex  drives  the  tail  force  in  the  opposite 
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COMPARISON  BETWEEN  PREOICTEO  AND  EXPERIMENTAL 

vertical  fin  normal  force 


direction  is  aetched  within  a few  degrees.  From  this  evidence,  it  is  con- 
cluded that  the  growing  vortices  tend  to  dominate  the  toll  ocrudynanic 
characteristics.  The  second  comparison  Is  shown  In  Figures  12a  and  b.  Here, 
predictions  of  Isolated  body  side  force  and  yawing  nouents  are  compared 
against  test  data.  Again,  the  effects  of  considering  all  of  the  wake 
vortices  separately  from  the  growing  vortices  are  considered.  To  make  the 
distinction  clearer,  no  attenuation  of  non-growing  vortex  effects  has  been 
included.  It  trill  be  seen  that  use  of  all  the  vortices  produces  divergent 
results  which  have  none  of  the  oscillatory  character  typical  of  asymmetric 
effects  and  evidenced,  by  the  data.  This  is  because  the  net  wake  circulation 
reshtii'S  unchanged  In  sign,  regardless  of  the  number  of  vortices  present. 

Use  of  ths  growing  vortices  only,  on  the  other  hand,  yields  quite  good 
matching  of  force  and  moo^ent  magnitudes,  as  well  as  angles  of  onset  and 
new  vortex  appearance,  at  least  until  several  vortices  are  present.  While 


matching  Is  not  exact  at  the  higher  angles  It  Is  clear  chat  use  of  Che 
program  will  produce  satisfactory  preliminary  design  level  estimates  for 
isolated  bodies. 

Lastly,  comparisons  are  shown  between  predictions  and  data  for  the 
10  caliber  body  with  a variety  of  tails,.  Based  on  the  results  for  tails 
and  bodies,  these  comparisons  contain  the  effects  of  growing  vortices  only. 

Figures  13  a,  b,  and  c show  side  force,  yawing  moment  and  rolling  moment  cor- 
parlsonr  respectively.  Prediction  accuracy  Is  generally  satisfactory.  In 
most  cases  the  tugnitudes  of  the  quantities  are  predicted  quite  closely. 

On  the  ocher  hand,  the  onset  of  asymmetry . and  the  appearance  of  new  vortices 
are  not  always  so  accurately  -eproduced. 

The  results  of  Figures  11-) 3 Indicate  that  the  vortices  growing  In  Che 
vicinity  of  the  body  dominate  the  Induced  effects.  This  appears  Intuitively' 
correct,  particularly  In  the  'context  of  the  two-dimenalonal  analogy.  There, 
Che  vortices  closest  to  the  body  tiould  be  expected  to  produce  the  pressure 
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rolling  mom 


•nd  velocity  distributions  which  generste  the  forces  end  moments , with 

i 

the  remainder  of  the  vortices  having  little  effect  as  they  pass  downstream. 

It  is  not  unexpected  then,  that  the  three-dimensional  case  should  give 
similar  results.  However,  although  the  street  vortices  have  been  shown  to  have 
relatively  little  Impact,  the  experimental  Information  on  their  strengths  and  | 

shedding  frequencies  remains  of  central  Importance  since  it  defines  the  i 

bounds  of  growing  vortex  strength  and  shedding  frequency. 

In  general,  the  model  performs  quite  reasonably.  In  view  of  the  changes 
in  force  and  moment  data  rsgnltudes  and  signs  which  can  be  obtained  by  rotation 
of  test  aK>dels,  some  degree  of  fortuitousness  might  be  assigned  to 

the  results  shown.  On  the  other  hand,,  it  is  unlikely  that  any  measured 
forces  and  moments  will  be  significantly  greater  than  those  calculated  since 
the  model  contains  all  the  essential  elements  of  the  vortex  flowfleld,  both 
in  magnitude  and  locations.  In  addition,  it  is  felt  that  the  means  of  calculating  | 
flow  field  effects  are  adequately  founded  in  theory.  I 
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6.0  PROGRAM  USERS  MANUAL 


« 
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6.1  Input  and  Output  Formats 
Input  Fonuit : 

The  Input  requirements  for  this  program  are  relatively  simple.  The 
input  can  be  separated  into  9 basic  groups.  The  following  gives  Information 
concerning  each  grouping. 

Group  I - IFlov  field  data 

a.  Format:  4F10.0 

b.  Variables: 

VINF  - Free  stream  velocity  in  ft /sec 
FSMN  - Free  stream  Mach  number 


RHOINF  - Free  stream  density  in 


lb  sec 
ft« 


> Maximum  body  cross  sectional  diameter  in  ft. 

.2 


and  - Kinematic  viscosity  in  ft  /sec 
Croup  II  - Body  Data 

a.  Format:  6F10.0.  213 

b.  Variables: 

D 

SREF  - Body  reference  area  in  ft* 

NOSEL  - Nose  length  in  ft. 

BODYL  - Body  length  In  ft. 

XMC  - Moment  reference  measured  from  the  nose  in  ft.  (absolute  value) 
DELTA  - Nose  half  angle  in  deg.  Delta  - tan~^ 

noss  isngcn 

NS  - Number  of  segments  into  which  the  flnless  portion 

of  the  body  is  to  be  divided.  (Maxlmum_ of  100  segments) 

IDCONF  - Configuration  Type  (0  • Body,  1 ■ Body  + Tall) 

Croup  III  - Fin  Data  (Input  only  if  IDCONF  > 1) 

a.  Format:  (7F10.0) 

b.  Variables: 
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XF  - Axial  distance  In  feet  from  nose  to  physical  fin  root  chord 

leading  edge.  Measured  negative  from  nose  aft.  (See  Figure  14) 

YF  - Distance  In  yaw  plane  from  body  centerline  to  fin  root 
chord  leading  edge.  Measured  In  feet.  (See  Figure  14) 

ZF  - Distance  In  pitch  plane  from  body  centerline  to  fin  root 
chord  leading  edge.  Measured  In  feet.  (See  Figure  14) 

j 

XFMAX  - Axial  distance  In  feet  from  nose  to  fin  tip  chord  leading 
edge.  Measured  negative  from  nose  aft.  (See  Figure  14} 

YFMAX  *•  Distance  In  yaw  plane  from  body  centerline  to  fin  tip 

chord  leading  edge.  Measured  In  feet.  (See  Figure  14) 

ZFMAX  - Distance  In  pitch  plane  from  body  centerline  to  fin  tip 

chord  leading  edge.  Measured  In  feet.  (See  Figure  14) 

2 

SREFT  - Tall  single  panel  reference  area  In  ft. 

The  above  Masurements  apply  to  a single  lee-side  fin  in  the  pitch  plane 
of  a non-rolled  missile. 

Group  IV  - Indicators 
a.  Format:  (313) 

. Variables 

NAOA  - Number  of  angles  of  attack  to  be  considered.  As  mary  as 

30  angles  of  attack  may  be  entered  per  run. 

KTYPE  - Nose  type  Indicator 

NTYPE  Nose  Shape  I 

1 Cone  \ j 

2 Tangent  Ogive  I 

NLAM  - Number  of  roll  angles  to  be  considered.  As  many  as  30  f 

roll  angles  may  be  entered  per  run.  j 

Group  V - Angles  of  Attack  | 

a.  Format:  (8F10.0)  j 

b.  Variables; 

AOAD  - Angles  of  attack  in  degrees  (See  Figure  14) 

Group  VI  - Roll  Angles 
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a.  Format : (8F10.0) 

b.  Variables: 


RA*  ■ Roll  angles  In  degrees  (See  Figure  14) 


*Locates  fin  positions  relative  to  leeside  vertical, 

RA  - 0.0 

Group  VII'  - Radial  Vortex  Limit 

a.  Format:  (FIO.O) 

b.  Variables: 

CAMLIM  - Radius  at  which  the  influehce  of  a vortex  on  the 
body  goes  to  zero,  (diameters) 

Group  VIII  - Hose  Fineness  and  Bluntness  Ratio  Option 

a.  Format:  (11, FIO.O) 


b.  Variables 


lOPTl  - Option  indicator 


(0  ■ Do  not  use  option,  1 • use  option) 

BRN  - Bluntness  ratio  in  percent  (s§fHSrfsif.  * 


Group  IX:  Run  Configuration  Indicator 

a.  Format:  (II) 

b.  Variables: 

IRON  - Run  indicator  (1  - another  run ' follows, 
0 ■ terminate  computations) 


Output  Format : 


Output  format  will  be  determined  by  the  nature  of  the  configuration  being 
analyzed.  Output  for  isolated  bodies  consist  of  angles  of  attack,  side  force 
coefficients  (CY)  and  vawing  moment  coefficients  (CETAMC)  about  the  moment 
center  input  by  the  user.  If  the  fineness  and  bluntness  ratio  scaling  option 
is  selected,  output  will  be  scaled  and  repeated.  Output  for  body  plus  tail 
configurations  consist  of  angles  of  attack,  induced  side  force  coefficient 
(CY),  yav'lng  moment  coefficient  (CFTAMr')  about  the  user  input  oKnnent  reference 
center  and  rolling  moments.  Additionally,  for  a fin  located  in  the  leeside 
vertical  plane  (RA  • 0.0),  the  force  normal  to  the  fin  is  also  printed  out. 
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6.2  Pro.^raa  Flow  Charts 

Two  prograa  flow  charts  are  presented.  The  first  (Chart  1)  is  a 
ganerallaad  flow  chart  designed  to  provide  Che  user  a basic  wap  of  pro> 
graa  functions.  The  second  flow  chart  (Chart  2)  is  designed  to  give  the 
user  requiring  a tiorking  knowledge  of  prograa  functions  a aore  detailed 
breakdoiqn  of  prograa  logic. 

Descriptive  stateacnts  are  also  contained  in  the  prograa  listing  in 
order  to  facilitate  tracing  the  steps  thn  gh  Che  prograa. 
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Chart  2.  Oatallad  Prcgraa  Flow  Chart 


A Raad  prograa  Inputa 


B Initlallta  calculatlona  for  a given  angle  of  attack. 


• .—1....  II—  I I .HI 

I Initlallxe  Indlcatora.  zero  oiit  atorage  and  calculata  axial  and 
croaaflow  coaponenta  of  flow. 

Locate  Thonaon  vortex  atrength,  aeparatlon  angle  and  aeparatlon 
point  paranetert  in  data  tablea. 

^T—  _•* 

I Calculate  potential  flow  generated  vortex  atrength  paraaeter.  | 

Calculate  vortex  atrength  acaling  factor  for  super-critical 
Reynolds  nunbera. 

~ i * 

Define  Scrouhal  nuaiber  as  a function  of  crossflow  Reynolds  nui^r 
and  configuration  type. 

^ i ■ 

Define  nose  shoulder  location  in  fon  of  separation  point  paraaeter. 
* ~ 
Calculate  Thoason  vortex  trajectory  angle  using  separatioa  angle 
paraaeter. 

— 

Establish  region  over  which  base  effects  will  influence  vortex 
growth. 

\ Calculate  vortex  clrcuaferential  separation  angle.  | 

I Calculate  vertical  distance  between  street  vortices  of  like  sign.  | 


Calculate  the  lateral  distance  separating  vortices  in  a stable 
Von  Karman  vortex  street. 


Scale  street  vortex  strength  paraaeter  for  super-critical 
Reynolds  nuabers. 

I Calculate  street  vortex  strength  using  strength  parsaeter.l 
I Calculate  street  vortex  growth  rate^^ 

Deteraine  effective  body  length  under  Influence  of  vortices. I 


For  body  alone  configurations  effective  body 
length  equals  total  body  length. 

* — 

For  body  plus  tail  configurations  effective  body  length  la 
the  axial  distance  to  the  aost  forward  point  bn  the  fin 

laadlng  edge. 


. 


Dlvld*  th«  affactlv*  body  langth  Into  tha  apaclflad  auabar  of 
aagaanta  or  coapucatlon  Intarvals. 


If  tha  affactlva  body  langth  is  aufflclcntly  long  to  ahad  aora 
vorticaa  than  tha  Thomson  data  tablas  Indlcata,  add  additional 
straat  vortleasat  ragularly  spaced  straat  separation  intervals* 


Using  vortex  separation  point  pameters,  calculata  axial 
separation  pointa. 


Exaailne  each  body  segnent  to  datemlne  if  Thonson  axial  separatloo 
points  fall  within  the  segment . 

" 'I  ■ " p 

Alter  locating  Thoason's  separation  within  a segnent. 
deternlng  if  separation  point  Is  on  Che  nose  or  cylindrical 
portion  of  body. 


I'sing  local  body  diaawcer,  calculate  coordinates  of  point 
where  Thoason's  vortex  trajectory  penetrates  surface. 


I Redefine  vortex  axial  vortex  separation  points  as  chose  points  sham 
I Thomson's  vortex  trajectories  penetrate  the  body  surface. 


CMrt  2.  (cotttimM^) 


EnaBln*  •«ch  t«ta*nt  to  dctoralmr  th«  nuabor  of  vorttea* 
on  or  prior  to  It. 


Dotoralno  tha  nuabcr  of  ahod  vortlcaa  paaatnc  ovar  Cba 
aaitaant  and  thalr  affect  on  tha  aagaant. 


|Shad  vortlcaa  pasalng  ovar  tha  ant  Ira  aaiMut.  I 


^ Ara  thay^»«i,. 
atraat  vortlcaa?* 


Vortaa  traj'actorlat 
daf load  by  Thoaaon 
saparatlon  angla. 


Dataralna  if  tha  vortlcaa  ara 
tranaltlonlng  froa  tb^  Foppl 
growth  trajactory  to  tha 
Thoaaon  trajactory  or  If  tha 
tranaltloo  la  eoaplata. 


' Circulation  about  aagnant  daflnad 
ae  coaponant  of  vortaa  circulation 
raaclvad  patallal  to  body  longltudl* 
aal  aala. 


Calculate  contribution  to  aagaant  aide  force  and 
yawing  aoacnt  accounting  for  the  decay  Ir  circu- 
lation with  radial  dlatanca  aaparatlng  too  vortaa 
■ and  aagaant. 


If  any  vortlcaa  are  ahed  within  tha  beunda  of  the  aigaint . 
dataralna  the  affect  of  tha  ahad  vortlcaa  on  tha  portion 
of  tha  aagaant  over  which  thay  paaa  taking  Into  accouat  as 
abova,  tha  trajactory  they  arc  on  and  the  decay  la  clrca- 
latlon  with  radial  aaparatlon  dlatanca. 


'V' 


Chart  2.  (ce«tlaiM«> 


Chairt  2.  (eontlaiM4>' 
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Chart  2.  (cMtiBiMC) 


U eha  jrawlm 

atraac  rarta*!^^ 


OaCanilna  If  tha  aagaant  la  oa  tha 
noaa  or  not  and  eaXculata  tha  ver- 
tax  growth  accordingly. 


Arc  tha  growing  wortleaa 
Initial  atraat  vertlcaaT 


Calculata  worta*  growth  aecotmtlng  tor 
any  baaa  affacta  and  tha  aodlflad  growth 
rata  raaultlng  froo  tha  irregular  apae- 
Ing  batwaan  aoaa  vortleaa  and  tha  Initial 
two  atraat  vortleaa. 


Calculata  vortaa  growth  ualng  tha 
noailnal  atraat  growth  rata  and  ac- 
counting for  any  baaa  affacta  on 
growth  rata. 


Calculata  contribution  to  aagnant  alda  foreo 
and  yawing  aonant  accounting  for  tha  dacay  In 
circulation  with  radial dlatanca  aaparatlng 
tha  growing  vortaa  cora  and  aagnant. 


If  no  vortleoa  ara  ahad  fron  aagnant  of  intaraat, 
calculata  contribution  to  aagnant  alda  forca  and 
yawing  nonant  aa  a raault  of  two  growing  vortleaa 
paaalng  ovar  tha  aagnant. 
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Claar  totAl  alda  forea  and  yawlnt  aoaaat  arraya  kafora  avaawtioa  for 
aagaant  of  Intaraat  baalna. 


Zaro  out  contrlbutlona  to  aid#  forca  and  yawing 
aoaiant  froa  al^  vortlcaa.  Thia  portion  of  pro- 
graa  aaoda  to  bo  raanvad  If  thaaa  contrlbutloaa 
ara  to  ba  eonaldorad. 


Sub  contrlbutlona  to  aagaant  alda  forca 
and  vawlna  aonant. 


Calculata  aagaant  aoaant  ara. 


Hova  to  neat  body  ae 


Sub  alda  forca  and  yawing  aonanta  for  all  body  aagaaata 


la  thla  a body  plua 
tall  eonf IguratlooT^ 


■ Call  tha  oubroutina  r'  calculating 

wortaa  affocta  on  flna  (riMCAL)  and  calcnlata 
I fin  noraal  forca.  Induced  rolling  noapat  and 
[ ^ contribution  to  alda  forca. 


Sub  laolatad  body  alda  forcaa  and  yawing 
aowanta  to  fin  contrlbutlona. 


[1 Calculata  alda  forca  and  yawing  noBant  ceaff Iclanta. 


I J Co  to  neat  angla  of  attack  and  rapaat  above  calculatlona.l 
K If  body  plua  call  configuration,  print  output.  I 

t * — 

K Tf  Indicated  for  laolatad  body,  Bodlfy  output  for  neaa  flaanaaa  ratio 
I and  bluntnaaa  af facta  and  print  output. 

[ L C.iack  to  dataralna  If  another  configuration  la  to  ba  eonaldorad.  I 
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C SUM  SIDC  FO-JC"  AND 
Wo  ll2  <al«^S 
TSF«TSF*SP(<) 

TYM:>hx  frMOH^Y  lOm  o 
3l?  CONTiNUF 

C GO  TO  '5S  IF  ISOl.AfE 

IF{1JC^njF,E3.'i)'?U  I 
JMAXi>JAAX 
DO  it>0  <«1*JMaX 
lF(AbS(fiSFPikj ) .l£‘. 
JMAXlxK-l 
GO  TO  is? 

COM  I 1N.J£ 

NFF1=<) 

JMAXsJMaXI 
FNFl tl)*0,0 
S»n.H{  l)ai»«0 
OaL-  FInCaL  I*»  CA-l-f 
CAU-  FImCaL(J'<Ax.  jo 
lGSE»ni)«SS£P7(T»tH 
?RNF.hYCP»QO,,M) 

IFdSTOP.EO.l)  ISFI® 
C IF  IST(J-»=1  A v/nor£.« 

C OISCUNTiNuE  CoLCU.u 

1F(ISTIIO,E0.1)GO  T) 
SPO,«n»*SHOL'»(T>*-T 
C SUM  SIOE  FO^C”  aNJ 

TSFaT  jF*W'jFy 

TYMJMeTYHOM* ( ( AnS(E 

*■  IF^.AM,Ew,0,0^N^^  Is 

IFr.AM.FQ.O.lDFMF  I ( 
M0RAsn0uA»1 

C Go  TO  S«  WHEN  CONT^ 

IF( yORA.GT. YL*M)&D 
I.AM3KA  ( vjO‘>A) 
LAmJ»LA4 

LA'AsLAM/Sr.?**^/? 

C CALCWLATE  CJOsOTNAT 

AFPsXF 

VFRs-l,*(0/?. » *Sl Y( 
ZFns«0/4.»*Cn=:<L'''’> 
GO  TO  JSl 

C CAL-OlatE  F:?fcE  sT-»E 

S8  0*,3*PH0I'JF*|./ImF** 
C CaLCOLATE  SlOr  ru^C 


IB- 


■iHTY  PRACIIC^^ 

, , 

OnOOYAAO 

) SI  )£  FOHCc  ^M■T  YAMTur  Om  Ea^H  Sc-GmE^T 

>*sFPSVt  jT»*srs‘jTS(jT>*SFiJUAS(  iTi*SF3BU(  |T) 

{ jT  ) ♦F  MtonSv  ( IT  ) ♦cMfjR.iT  j JT  ) aEmB  ia<  ( JT) 


iNIEN  'TF  P*Tt«S  iMe 
5(sF«K)  » 

ANji  REPEAT  CAl  C'Ji  ATT'AMc 


YAWING  MOMENT  “O'-TR I »«IT  1 MNS  FFnM  ALL  SEGmEvT^ 

A-. 


) iOljY 
D bS 


\\ 


ARsIEHOUTL)  T-i  350 


3Sb 

352 


351 


% 


^ODYL)  ♦ TAbSTxr  4AY-XF )/■>.)  )»HNFy) 

1 

I ) *4  -IF 

iBJTiONs  fwum  ali  firs  Have  reen  calculated 
rj  is 


->  iF  FlJ  ROOT  C-''|WN  lEAOIi^s  E0';E 
,A'^I 


A4  DYNAMIL  pJC^*;  H>C 

?.) 

■ Ago  Taring  -"iMrir  '•ncFFlClENTS 
69 


0o007730 

On0077An 

0n0C7750 

0(1007760 

On00777o 

On00778o 

0o00779n 

0n007P0n 

OaO0751(1 

On00782o 

0(i007830 

0ii0078Af> 

0n007a5n 

0(i00786o 


-i-t ' 


On007ASn 

0A007A60 

0n0O7A70 

0OO0748O 


0(1007490 
0«00750n 

00007510 
0o007520 
00007530 

0n007540 
0O007550 
0000756(1 
Oo00757o 

0n007580 

0o007590 
0n007600 

0o00762n 
0o007630 
0O007640 

0o00765o 
0fl007660 

Y Jut  VORIEA/PTN  INT•9^CT ION  CilrULATlOMA 
■lb<l(l)«XFKtVF4,7F4,xPM/(X,xFt7FiAv»£F*LAM»  Oo00767o 

«t  tol  » J»  vInF  a «S9e’FT»9HolMF  tG5EB(  i )»HnFy»tSTOs*  0(io0768o 

0O007690 

(I)sl  O0OO77OO 

OAsStS  aITHT-  1/4  SOOv  OAOIJS  or  FJm* 

IlDNStjU  To  -JC^XT  ANRlF  of  attack 
oi  0o0077ln 

.»LM  O0OO7720 

TARING  RiOMENT 


CY(I)«TSP/( 

CET»U)«TyMDM/(n*Si» 

tra'isfer  of  y»»iN3 
CET*HC(I)»« (CTTa<1) 
> continue 

IF(1DC0nF,E3..()5‘»  T 
PRjyT  body  bl  iS  T*I 
IF(yFFl,NE.O) wRiT- ( 
) FORRAT(3X*15HaNQL£ 
2AX*l2MNCtRMAL  •‘OrCE^ 
36X.»9H(0EUaEESt  »l*x» 
Ai6H(FlN  NO.  I (lB) ) 
IFtNFFl.Eo.O) -RlTE( 
) F0R'<AT(1X,15H4NGL£ 
r 6X»9h(')ESRECS>  *UX 
DO  313  I»1 tyAOA 
IFllSKlP(I)  .Ntr.o)  30 
IF<SFF1.NE,i)>  rfRlTEC 
9 FORNAT (7X,F&.?.6X«r 
IFCnFFI.EI,  J)  <<RlT- ( 
I FORMAT (7X,F6.Bt6XtP 

i continue 

60  TO  2a 

PRI  NT  ISUH  ATED  qOJy 
1 NRITE (6.2S) »M- 
5 Format (3X, 1 5MaN0L£ 
?6x«9H<0EGaEES) .lAX. 
00  26  lal.NAOt 
NRITE(6.27)A04D(I). 
7 F0RMAr{7X,F6,3,6X.e 
s continue 

IFdOPl.NE.DSO  TO 
IF  lOPT.I  scale  is 
0LUM'TnESS  RaTtO  EPF 
CALL'  OPTl  (FSM^||N0SE 
>IPASS.VA0<H 
IF(lPASS.EU.O»Gn  TO 
MRI  TE  (<).«nO) 

H format ( IX, IrfnalCK 5 
BRrTE(6,25)XH- 
DO  801  ralfMAOA 
RRlTE(6.27)a0a0(I». 
I CONTINUE 
A REAO  IHJN 

ChECX  to  see  tF  AMO 
IFdRUN.Eo.DsO  TO 
CALL'  exit 
End 


EF»0) 

MOMENT 

/CT  { I ) )-(XMC/nn*cy  iT» 

y 

.1  O'jTPUl 

OF  ATTAL<»AX,?HCy.6X.6NCFTAMC#7X,'9HR0LL  mOM,. 

^ t 

3H:vi»»F‘»,2«2MPT»«;x*7MIPT-L8)  .♦«§ 

) 

S.ln)XML 

OF  attack#  AX, ?-(Cv.6X.6mCETAMC»7A,9MN0LL  mOM., 
»3iCQ»#^5*2»?«PlT,SX.7M«FT-LB)  ) 

TO  313 

5#*)  AOAU(l)  .cy  (1»  .CETA**Cj1)  •SHOLm(I)  ,FnF1  (I> 
7.i,SX*F7.3*AX.F7.3#7X.FA.3) 

Stii ) AciAO(i)  ,ry(Tj  .cET^Hcd)  .SROi  m(i» 

7. J.5X,F7.3#A*,F7,3) 


OJTPUT 

OF  ATT4L<#AX,^^Cv.A<,XMCETAMC#/. 
3H:Ga»Fb,2,«iNET) 

:v { I ) »CfcTAMC ( T ) 

7. 4,5X.t>  7,3) 


OLATEO  BOOy  OllTB'iT  FOR  NOSE  FInEmESS  RaTtO  A> 

icrs 

Oi,  MRNf  xHC*ry  < 1 r *Cfta  f I ) tCETAwc  ( 1 > ♦ 


OPTION) 


:Y( I ) »CtTAMC( T) 


ThEw  RjN  is  TO  Ae-  M*-*? 

A] 


00007871) 

OnOQ788o 

Oo007P9f) 

0n007900 

0a007910 

0n007920 


On00797o 


OoOOROOO 

oooopoio 

OoOOROZo 

0oO0A03o 

OaOOBOAq 

OoOOrOSo 

00008060 

0n00s070 

onooposo 


ohoopiio 

onoosiZo 

OnoOAl3o 

OoOOAlAO 

oooos;5o 


00008160 

OnoO0i6i 

OoOOR170 

OOOOA180 

On0OA2Oo 

OnoOaZIO 

OnOOn?Zo 

On00323o 

On008Z7o 

OnOOSZSo 

OnOO029o 

0O00R300 


iHIS  PAM  IS  BEST  QUALITY  PRASIIUUJUI 
fBOM  CX>Py  FUfWISHED  10  MC 


no  onoooo 


Ail 

SumOUTiNf;  KY7(hO»  THiir*!  J'4-»««,Y.T«Ol  ♦OELTAi 


SUB^OUTT^E  Ail  CA-CJL^TES  CoOpaTMATr*^  OF  pOT^TS  AlOmG  V^HTF* 
eoorfTH  T«*J£CtOpY 


IF  UBGCaJ)  , 5T.  AqS(»  JJ  » f;o  TO  1(10 

Calculate  cyopOtnatib  oe  points  \.0m«;  «qowT‘4  tra^'fctOSy  of  a 

MOR'IAL  STPEiT  V()HIF< 

PI  aATA-^t  ( (^/P,)-(  ( VE,»*C05(TMFrali  ) »/( AAS(Xj)-AASr*JB?)  ) j 

CC«AB*(AHS«<J)  -X»/:05(PI) 

ACBAH  ■ Si)HT(  ( ( (M/?,) -(  {0/2,  > *Cn«;{T^FT»T  > > )*»2  > ♦ 

J { (A3S{*.))-APS(XjMi)  )••?  1)  , 

Imti/?,)  - (l*8^(XJ)*  X)*TAn(»1)) 

101  delta  « ATA>4(  (SqhT  (((  i/a.  » ♦SOPf  f c (-</'». »•••»•)♦(  (0/?,  >**2,  » I ) aAE, 
1-^  ((•»/?.»  ••2.) )-(  {J/’.)*5Is(THEt  »T)  1 »/ACPAo) 

Y»  ( (ACBaM-CEBaP)  •{  AYi(OtLTA»  )♦<  {1/’.»«r'*«<THETAr«  * 

SO  TO  102 

Calculate  CDo^oiNAtifi  »e  points  a.o*>r  rpo«..t*4  tmajcctOpy  of  a 
VOPTEX  FOPCiO  To  5E3a-<aTE  mT  Tmc  aa«;F  nc  tHF  900Y 

loo  Pi  « atan(  ( (H/2,»-(  (j/?.)»cos(  rMFTAtm/(APp(Aj)-AMp(xjA2jn 

CEBAM  »(A9S{A  |)-X) /:oj(BI» 

ACBAR  aSOPTC  ( r (h/2,)-(  (a/a.lACOPTTMcTAH  Jl|A»2  )♦ 
lt(A3S(Aj)-A3S(X.|M2»)«*?  )) 
z«  (H/2,»-((AaS(AJ)-X)aTA<{01»> 

GO  TO  lol 

102  BETJNN 
ENr» 


AnOOBBln 


0ft00A32(l  ‘ 

I 

j 

I 

OflOOP33n  I 
0a00a3Ao  I 
Oflo6«3So  ; 
0A00B36n  1 
OaOOaSTo  1 
0O0OB380  ! 
0oo08390  I 
OaOOaaOo  I 
OnOOBAlO  7 


0ft008A20 
OndOSAJO  1 
OfiOOBAAo  { 
OflOOBASO  ] 
OnOOAASO  i 
OftOOBATn  I 
OnOOHiSo  ] 
OnOOflABo  i 


•t 

I 


E3P> 


SUBROUTINE  £090(1  NT* I -O)  OnOOBSOQ 

• #♦•••••••••••••••• •••»•*••••••• ■ llw 

subroutine  £0,0  DETiRNlNES  A VOBtE*  TS  NuMHEbEo  rVEN  Or  nOO., 

This  OETERMINrS  T-»E  Sl3N  4PPLItr»  TO  THF  CTBrULATlriN 


InTJE^InT/s 

INP1»INT*1 

lNPl02«INOl/2 

IEB»i 

ir(INPlr)2.6T.TNT02»  IEJ«0 

000  ICObO 

EVEN  IE0»1 

RE  I jRn 

EHO 


OfiOOnSlO 

OnOOB524 

OO00B530 

OftOOBfRo 

OoOOrSSO 

0A00B560 

oaoorsto 

OaOOBSSO 

0A00B590  I 

I 


SUB-<UUTINE  pI\jCaL(  A4,  OnOOAAOo 

I USE®  I t'»SEPZ»H,E»PHI  »0»vlNFMtSWEFiPH'tlMF,RciP»HNFy»l«TOP,  0n00fl6l0 

?RNr»WYCP»POtM»  OnOOA62o 

REAli  LA‘4  OnOOA63o 

Dimension  gamsviI) . ssiPKi) tOSEP7( i » #Gsrp»i» ♦*<ionv(.' o) tzcioj • OftOO«6*o 
1 VMT { 10) » VT ( 10) » AOA05 « In ) * Vn^AP (?n) tcN tP>)  OnOOASSo 

dimension  YCP(2n)  0n00A660 

subroutine  fImCal  CAL-tJLAltS  IMF  ^0->'Ui  Fn«fE  InD  jCrU  On  FInS  Br 
shed  VOpTiC-S  IN  T^E  HOOY  i.tE  SIDE  M*<c 

JST4RT»,)MaX-1  OoOCRBTo 

Do  53  JsjmART*  |M4»  0«OOR68o 

shed  VOpTfX  StRcniTHi 

GAM*0AMS'/(  J)  OnOOAG^O 

CALCOLATlON  Or  cODRUNATEi  ON  Fin  LCAOt'iR  Eo6£  #rtrpp  CA.CulaTIOnS 

ARE  TO  PE  MADc 

X(1)»XFp  OnOOftTOo 

Yn)=YFp  OnOOPTlo 

2 ( 1 ) *ZFo  OnOOflTBo 

FLF.S=(XFMAX-Xr) /(Zf>tAX-7M  OnOOfitSO 

0ELZ*<^FNAX-2r) /lo.  OflOOBTAQ 

DO  51  I*U10  OnOOPTSo 

IFd.GT.DGO  TO  56  OnOOpTGO 

X < I ) =FLrs» ( DE .2/2, ) • X ( I ) OnOOPTTn 

Z ( 1 1 =Z ( I ) ♦ ( OE  2/2, ) »C JS(LAM)  OnOOATSo 

Yn)=Y(I)-(DElZ/2,)»SlN(LAM)  0 no  Of  TOO 

GO  TU  S7  OnOOpBOo 

56  K«l-»  OnOOBSlO 

xn)*FLES*UELT**(<)  OnOOpftEo 

Z(I)sZ(K)*i)iL7*r05(.»^)  OoOOB83o 

Y(l)=yi<)-U'LT*<;M(.AX)  OnOOfiBAO 

AXIAL  OISTANC-  FHD*»  0Jr>(T  ur  YUPTE*  «Eo'«»TTON  TO  PnlNI  On  riN 

57  0ELXaAHS(X(I),GsE3T( J) ) OnOOPBSo 

definition  DF  VOPTEX  .IaTEPAi.  OISP.Af'FMFNT 

IF(3AM,GT,0,0,AnO,  J,Li.2)  YVCL»-r,S-P7(.J)  OnOOgSGO 

IF  (3AM,Rr,0,0.  AmI),  j,r«T  .?)  YVCL*-H/a,  OnOOPBTn 

IF  (3AM,GTii),0,  ANU,  J.E  J.  jSTAPr  .OP.  j,  r -) . .|vi  a » ) y VCL  x-GSrPZ  ( j)  OnOOPBBO 

IF(G4M.LT,0,n,AN3, j,.E,2) tvLL=Oc-P’(J)  On008B9o 

IF(3AM,lT,O,0,AnU,  J.3r,2)  YVCLo-^/P,  OnOOBPOo 

IF  (3AM,Lr,0,0,ANO,  J,E3.  JSTAPI  ,0P.  J,FO,  J'*A*  ) YVCLsGaEdZ  (J>  0n00R9l0 

AI»  ( ABS  (2  ( I ) -oEl  '*T4'4(r  ) ) ) *CUS  «F)  0n00R920 

Bl=A3S(yYCL-Y(I) ) OnOO«93o 

ANr,.E  = ATAN  ( Al /H)  ) OnOOB94n 

a8AR=0ELX/C)S (E ) 0n00p95o 

OBF»ARs(APS(DE.X*Tf  N(F)-Z(I)  ) )*STf'('')  Oo00b960 

BBApsPOAP-UPHaP  OnOOfl97o 

ZAO  J=‘8BaH»3iN(t)  0n00898o 

Phi  1 =A-if'.Lv:  OnOOB99n 

PaoIAL  Tl‘;raN-E  FP')'  /Iprrx  cone  TO  POT'iT  ON  FIn  viF,SJH-0  m'OnG 
A LiNF,  =>E-'r-‘?:N-)ICO.A  / r.)  Ht  VOPTFX  -obf 
67  »I*  ( AHS  I Z ( I) -nEt,  X*T  AN  ( r ) ) ) ACOS  (Fl /S  t (OP1 1 ) OoOOpOOO 

P10*p  1/ {U/Z.  ) Of)O09OlO 

IF  yOPTE*  paSPES  AJTPI  I 1/4  body  PA-MiJB  Or  POINT  ON  FiN 
:/is;>jNfiN)E  c^Lc-'-A^n-js 

?r  f »in  .UE  . 0,2^  ) oO  T I ?1  0f)009020 

vEi.^CSti'  r«JU'E)  A ?■  =)INT  04  FIn  ^y  VOPTEy 

YT » • jap/ f ?. •p^ r •“ i ! 0oo0903o 

t,  )F  YrLDtITY 


73 


0f)009020 

0n009030 


VYl»'/ri*SlN(P-iIl  ) 
VZi«VTl*C^S(P^Il ) •CDStE) 
lFCZAaJ.I-T.2(T)  )GD  ro  t<t> 
ir(rVCL,LT,r(T) )6D  ro  M7 
VY1»-I.»vyl 
GO  ro  69 
A7  vri«-l.»vyi 
VZ1=-1.*VZI 
Go  10  69 


86  lF(yvCL.LT.r(T) ) VZ1=-1.*VZ1 

69  VYNi»vr 1 *cu; (lAh)  , 

VZn1»VZ  j*Sl  Yt  >9) 

N0P9AL  CO'^P.INENT  JF  Vil.uCllY  lNOi)';e'^  OM  FtN  BY  YOBTrK 

VMjeVYNl*/ZMl 

AZ«AaSt7(  I)-DpLX«IA'Ji(*')  ) 

B2»4dS(YVCL-Y(l) ) 

/^NGLE2»ATaN  (A?/82» 

ZaDJ2»0ELi>  »TaM'P) 


PHI  2 

70  9?a  ( AHS  (Z  ( I ) -3Fi.*»T  a V ( f ) ) ) /SiM  (P-tt 

IF  ( yVCL.GT.Y  ( I ) . and.  za  )J2.6T  .Z  I T ) » 5/^  TO  «« 

IF  ( YVCL.LT.  Y ( T ) . AND.  Za  .Z  ( T I ) So  TO  Ao 

IF  ( YVCL.PT  . r < T ) . A \'0.  z ao  j<>.lT.  7 I T 1 ) So  TO  Bo 
IF ( ^ YCL.UT . Y ( T ) , A NO. / a )j2.t r . Z ( T ) > So  TO  Oi 
OFFlNJTlJ'i  Jf  EoUI  VAL  i'jl  ■ VOP  1 F A STPcA'fiT-l 
86  GAm2  = V\1*(2.*3MI«-IZ)  / ( O 'J6  ( H-i i ?)  *<:  j n f L * h) -aI N (Phi  2)  arUS  ( .AM)  i 


GO  TO  70  , , 

89  GAM2»-l.*vNl*(2.*3Hl*NP)/(SlH(Pm3).rns(LAM)  ♦00$  (bHt2)  •SiMIi  AM)  i 


SO  70  73 

90  S4M?  = VH1*  (Z.*3HT*9r')  / t S I N ( P-t  1 £> ) *0  OS  ( I.  AM » ♦rOS  < PH  I 2 ) •«  I M ( .AM)  J 
30  tU  73 

91  SAmZ»V^J1*(Z,*ami*9Z)  /(SlN(HHl?)*C0S(LAM»-rO$(PMl2)««IM(.AM)') 
LOCaTiON  OF  I-ISE  vOalcx 

73  F«saBT  (((  465  YVrt »)  •*Z.  )♦  (UE  LA  *7*  N(c  3 ) ••?.) 

QI« { (0/3. I ••?. ) /F 

THE  r A*  AT  4^1  ( (A4S(YvC-»  l/(UFLX*rAM{r)>) 

IF(3AH?.Gf.n.1)  rH;r4*-l.*rHi:-lA 
Y I»  31  ( T HF  TA  ) 

ZI«3I •€ OS ( T iE T a ) 

PlaSUHT(((A3S(Y(I)-Yl))#*4,)*<Ao<;(7(I)-7T»)**2.) 

GAMl*-! .*54^2 

CAlC'-IUATE  YELOCiTr  I'jJijCFO  af.PoTNT  ON  CIM  qy  I'^AOE  VOHtE* 

YTI rSANT/ (2.*3HT*9I  ) 

••21  sAHS  (Z(l)-xl) 

821 sAHS( r r-Y ( T ) ) 

ANOlE I»4TaN ( A3l /8?I ) 

Ph12Is4nOlEI, 

VYI»VTI*SIM(P-'I?1) 

VZl«vTI*C:)S(P-ir3l) 
iFtZl.LT.Zd)  )G0  1'^  yZ 
IF(ri,Lr.Y(i);6n  TO  93 
VYT*-l.*VvI 
60  TO  76 

93  VYI«“1.»YYI 

GO  TO  7y, 

92  IF (Yl.LT.v . n : VZls-j . *vZI 
76  ¥YNI*VYI*C03(..AM) 


OnOOqOAO 

00009050 

O0OO9O60 

00009070 

00009060 

0n00909o 

OnOOqlOO 

O0OO9IIO 

0n009l20 

Oo009130 

O0OO9I6O 

00009150 

00009160 

00009170 

00009160 

0fl009l90 

0ft009200 

0ft0092l0 

00009220 

0ft009230 

0OO09240 

00009250 

0o00926n 

0o00927o 

0o00o280 

00009290 

O0OO93O0 

O0OO93IO 

O0OO9320 

00009330 

0o0093AO 

00009350 

OnQ0936o 

O0OO9370 

O0OO9380 

Oft00939o 

O0OO9AO0 

O0OO9AIO 

OftOOqA2o 

Oo009*30 

0O009AA0 

0o009A50 

O6OO9460 

Oo009A7o 

On009ASo 

Oo009A9o 

O0OO95OO 

O0OO95IO 

00009520 

Oo00953o 

O0OO9540 

Oo009550 

O0OO9560 

00009570 


^i 


uuuuu  u uu  u 


PI 


W7HI*VZT»SIV( ,*-) 

>iO0<4AL  Cv)'4PDNsimT  DF  ViL'JCIIr  Ii^Ou-Cr,  *T  P^I^f  0«<  riv  it  I^a^E 
voprE* 

VMlaVYlMl*V/^I 

RCaSUPrcAPSlYdi  )••?.»  ^^S(Z(  I) 

CALCULAflON  Op  CO'<‘'ON£nT  1C  »PL«»'*lTy  INO  iCn)  AT  aOpiT 

ON  ‘flN  PY  C£NTH»L  YO^iIEX 

VNC»GA'«?/<^.*3Hia^C) 

SUM*4ArinN  or  ^aOPH«l.  viLOCllY  COmoimcnTS  p«P  all  SrbjCNT? 

VNT ( I ) *YNI*VNl  ♦V^C 
51  CONTINUF 
VT(.MaO.O 
00  <1:  lal.ia 
YT(  WaYT*  J»  *V>iT  I I > 

5?  COnyI'^'f 

00  i-i  1x1.10 

effective  avG.F  0“  •rfaCK  iNUjtCi)  4r  POTN»  on  FiN 
AO  A 3i(  1 » *AT4N  r A4S  ( W 'IT  ( n > / V li'(F«» 

IF  t VNfU  » .Of  4 HT>(  I J «-l  , t » 

A3  CONflNUf 
5Y43Yxa.O 
SAI)Y«n,ft 

00  Ixl.lO 

IF(  I.EO.l ) Y:PTxn/^,*OiLZ/i-. 

IF(t.or.l»Y:PT»YC»IO£Li! 

SYA3YxAr)A'>S(I»ArCaf  •o:lZ*<T40y 
SAOY*AUAOS(!»*OrLZ«5A3T 
NA  CONflNOi 

F|n  SpA'iNISE  tE''Ti:i  0»i  apSbSOpt 

YCPI  J>x^YAOr/^A-ty 

4VF-4AOt:  VFl5CTTv  4I  Fl^t  I P4-NtN0  Ff*oF 

vnpa>4  (3)  ^AtlS(  </t  ( J)  / t I.  n 
AlN'»Ar4M(  V Jl  / Vl  tFA) 

YFL-«“S'4PTMv1''F4*»?.I  •(VNHA-^(  I 

AP»?.»  (ZF '•A  4-rFj  JEF 

CnAsPHI *44/^. 

3»0.»*’<‘<OtNP»  ( VFL  y*»?.  ) 

NOP  »«L  F0-»CE  OY  O V'>TF« 

FN(  J)  xCxAaAIN'jxoaSJ-F 
AlMl*Al*»0*5f.5<» 

IF( ¥ r ( J) .lT.O.OiF Ml J» x-i  .•F  j(j) 

Lon’  AND  CA.CiLuT*  mOmhaL  P0*<CE  (mO  -CFO  By  »LL  OT^Ea  V0?Tlrr5  tv 
04^1 

53  C0NII"*UE 
Ri^FxO.a 

PEBJLjAvT  FIN  NOH'irt.i  -O^Cl 
Do  5A  IxJSTANTt.jHAii 
RnFxPNF.Fm ( I ) 

5A  Conti MUE  ^ 

RNFr*HNF*CO^ ( Am)  -y 

rnf/»«nf*<;im( 'am)  ^ / 

PMXO.O  X 

calculate  HOL.J'JG  moment  ^ ^ 

00  95  IsJSTaMT.  i.max 
PM«TM*YcP(1)*pN(I) 

8S  continue 

ROLNB-l.ApNrA^YcP 


OaQOQSSO 


0(t00959n 

OA009600 


00009610 

0a009620 

0a00963o 

OaQOoOAO 

0o00965o 

O0OOQ66O 

00009670 

O0OO9680 

Oo00969n 

OoOOqTOo 

OflOOoTlO 

OflOOoTZn 

0o009730 

O06O9740 

60OO9750 

0ft009760 

O0OO9770 

00009780 

O0OO9790 

O0OO9OO0 

O0OO98I0 

Oo009020 

00009030 

00009840 

Oo00985o 

OoOOoOOfl 

0o00997n 

0OO09880 

0oO09090 


O0OO99O0 

0A009910 

0nO0992o 

0o00993n 

O0OO9940 

0o00995o 

O0OO9960 

0o00997o 

0n00998o 

00009990 

OoOlOOOO 

OoOlOOlO 

OoOlOOZo 


1 -.4  ••  ^ 'Y*.' 

■'w.r. 


SO  ro 
isr3**-i 

23  «eTJ«N 
iHO 


OoOtOOSo 
onoioo^o 
OnOlOOSo 
OaO 10060 


U»Tl 


SumOoTlNF  * >(»U»tlBN,A  tr  , C v .Tf  r « ,CP  T A <r  # T P vA3*» 

OI-^ENSMN  F^K  (h,<.)  , TM»(;H(^)  , IKf<  (4)  ,HRt'2(6»';)  , 

?flRK3(^f  b)  tB-»n*(».t3)»cr(i).CiTAi  H .cfi  AM-'n  > 

SUP^OUTINE  DPtI  SCft.Ea  ISOLaTEU  .oOv  SI-)E  F'^KCE  A'iD  rAKiKG 
MOMENT  FOP  FImEsEsS  PATIO  AMO  dMiv|TMF<;i;  p»TlO  tFFccrS 


OaOIOOTo 

OnOlOOSo 

0fl0l0090 


OAT*  TmaCm/0.  0.  ; 

DATA  TFM/?,0#3.0»3.7b9,A,0/ 

DATA  FM</i).  267  • n.  Si  1,081  * 1 . 037 . 0 , *0  3,  " , aBq, 

? 'i,933»0i  VA9*  I ,339»  1 , Adg*  A,17.T«'»,  t po» 

T i.nnn,i.oo'),  1.000  «i,ouo.i.nni)«i.Ano« 

A l,01*,n.<S7»0,M06f(',7l4.l  » .380/ 

Data  THj/p.oio.stin.iit/o.onO.o/ 

OAT*  MMH2/i,o->P.l.»no.i. 000*  1, UA«, 1.1.00,1.000, 

? l,U0.t,  1.000*  1,000*  1,00  1,1. 000,1.000, 

3 O,l-)0.n,??A*o,23M,O.  J41,0. ???,?. 7S0, 

A 3,24l,/,l 3A*?, lPA,3./S7,?.H1S,?.B7b* 

5 1 .Ood.O, 7 7A, 1 , 1 1 *, 1.341.0.2??.?.2aO/ 

DATA  MP^J/UOoO.l.UOn*!  .000*1, Uoo, 1.000,1.000* 

? I,0o0,l.i>)<)*1.0>l0*1.0oo,1.000*1.00p* 

3 0,703*0,A4S*O,V*O,<f.A4.>.0.r?O,l  .6*.  7* 

A 0,7  36,0,5??* 1 ,ASA*2.0oo,?.q7o,l  .778, 

A 0. 7AA.(t.  A 41  *1,  t.  I ;uno,p.  778/ 

DATA  PW<A/1 . Ono, 1, 01 .» *1. on 0*1,00, 1.1. lion. i.OoO* 

? l,OoO,l.oo<J*l,Oiiu»l,Oon*l.(>no,1.0oO* 

3 A77,0,S7S*n. 600,  J.  *40, 1 .083,0. S0<?* 

A 0,A77,.|.S»S*  ..61  J, 2, 330,0.213,0.176* 

A 0,?Pb*ii.‘,A|  *0, 700*3,640.  I .41  3,0.376/ 

IF  (FSM'*,I.  f ,'l,  A,<1M,C  ; *1.fiT.  1 , 1 K’6  rj  1 
FRU  = AM/I1 

IF  (Fkn.i.  T.2,o,0o,  • 8 ■ , , r , 4 .0  I ,V)  71  i 

1F{86m,gt,3o,'»)';u  id  1 

IPaAS»1 

DO  ? 1*1 *3 

J*  I ♦ 1 

IF  ( ',t  .TFP  (I».ftN),7|<i,l,t,1FHf  H)r.n  T6  •» 

2 COVllNiiF 

3 IFP»I 
IFP1*J 

on  4 131*6 
J*I»1 

IF  ( , TM«C8  ( 1 ) , A I.J.I  SM  4,  tt  . T ■,4>“M1  .11  ) rO  to  q 

A,  COM  TIM' IF 
5 IM*I 
1M1*J 

Ratio* ( tmacmi tm) -- s in j / ct^acm ( Imi .t  <achi 1 4i > i 
F I.fhk  I I I.  I=-h,  _,maT  T :)»  (F8<  ( IM.IFU,  -rPK  ( 1'*l  , IFP>  i ) 
F2=7H<  tiM,ls-W))«(^iriJ«(F3A(lM*fr3l,-F3<{|Mi,IFvil,>, 
FPK1»F1-(FI-F3)*( (T‘PIiF8)-fPnI /f TFQ(IF3).TFR(lFrti j j ) 
IF  HMM,lE,5,()i  6n  n ■* 

Do  f I*?*A 

J«T*1 

IF  1 5KM,3t  , f 8W  ( 1 1 . AN5i.  iRN.l.A,  f PM  ( ii  ) -.n  fn  * 

7 CONTlNUF 

8 1BR=1 
IBR1*J 

Fil»3HK3(lM,IBP)-HATl3*i(8PK3UN,IAR»-BR<3»lMl*lPF<i  i i 


.»  V 

*>  s 


OoOlOlOO 

Onoloilo 

onoloizo 

Onoloi3o 

OnOlOlAO 

OoOlOlSn 

Onolol6o 

0o0l0l7p 

OoolniBo 

OnOlOl9o 

0o0i020n 

0o0ln2l0 

0o0ln22o 

0o0l023o 

0o0l02*0 

0n0l0250 

0o0ln26o 

0o0l027o 

0o0lfl?8o 

0o0l029n 

0n0l030o 

0o0l03l0 

0o0l0320 

«00l033P 

OoOlOBAo 

0o0ln350 

0o0l036o 

OoOlOJTo 

0o0lft38n 

0o0l039o 

OoOlOAOo 

OflOlOAln 

OnolPA2n 

OnolPA3n 

OnOlOAAo 

OoOIoaSp 

OnOlOA6n 

OnOlOA7o 

OnOlftABn 

OoOIoaRo 
OoOlosOn 
. OoOlOSlo 
OnO I 052o 
OfiOl  053o 
OnOlOBAO 
OoOlOSSo 
0o0l056o 
0o0l057n 
0o0l058o 
0n0l0590 
0n0l060o 


OttTl 

(IqRi  ) ) ) 

ri»*^il-(Fl-F2)  »( (TiftlHR)— IHV)  /(TM^  C»hH)-to»(I8hi  » ) i 

ir(rRH,QE.3,0,AND.r?N,LE.*.0»SU  Tl  10 

82»aRK2(lMtl8:»l)-(R»Tl0«(8HK2(lM.T‘<'>n-4P.r2(l8l.I«Rl » ) » 
RATlOi»(2.0-FRN)/(2.0-1.0) 

00  TO  11 

10  6l»3RK4(m#lB?)-('<AriD*««RK4lIH.|lR)-RR<*jlMl»lPR)  » 1 

*i2»3RK*l  iMtlBpll-lRUlOPlBHK^dR^T-lol  J IMI  , iPPi  » ) ) 

RAT lOl ■(♦, 0-F PN) / 14.0 -3.0) 

11  ei*3i-<Sl-62)»((T3«(llR)-8HN)/(TqplTPB)-TpR|IBRl))) 
BRKlPOl-l'il-FnPM^TlUl 

60  ro  12 
9 BRRl"!.'’ 

12  DO  13  UUNAOA 
cy(i)«CY(i)»FP);i»3fl<i 
CCT*  n ) -CFT  A ( T » pF-IK  I • 3«<  1 

CETARC(I>«(  ICpTaII)  /Ct  (I)  )-(*'<C/nj  )*rv(M 

13  continue 
00  TO  U 

t IPASSpO 

14  RETJRN 

tNO 


OflOlOOlO 

0p0l0620 

oiiolooso 

06010090 

OoOlOOSo 

OnolOOOn 

OaOIOOTo 

OnolOOOo 

OnOlOOOo 

onoioroo 

onoiorlo 

Oa01072o 
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On0l07PO  i 
00010750  I 
0n0l076o  ’ 

OoOl077o  ; 
0n0l070o  I 
0a010790  • 

OaOIOBOo 
OoOlOBlo 
0a010820 
OoOlOOSo 
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6.4  SMple  Inputs/Outputs 


Isolated  Bodl«<i 
Input : 


VINF 

890.0 

ft /set- 

FSMN 

0.8 

RHOINF 

0.000642 

lb  aec^/ft^ 

ANl' 

O.OOOS75 

2, 

t i /sec 

D 

0.112 

ft 

SREF 

0.0768 

> 

ft 

NOSEt. 

0.936 

ft 

BODYL 

3.12 

ft 

XMC 

2.0 

ft 

■ \ 

DELTA 

9.45 

(IAk 

] 

NS 

50 

IDCONF 

0 

NAOA 

15 

i 

KTYPE 

2 

, 

NLAM 

1 

AOAD 

24.,  25., 

26. , 28. , 30. , 32. , 34. , 36. , 

i 

38..  40., 

42.,  44.,  46.,  48.,  50. 

RA 

0.0 

A 

gamlim 

100. 

lOPTl 

1 

• 

BRN 

0.0 

IRUII 

0 

Output : 

. 

ANCLE  OF  ATTACK 
(DECREES) 

CY 

CETAMC 

CC-  2,00FT 

24.00 

0.486 

1.209 

25.00 

0.342 

0.891 

26.00 

0.155 

0.537 

28.00 

-0.287 

-0.255 

20.00 

-0.791 

-0.909 

1 

32.00 

-0.705 

0.308 

34.00 

-0.241 

2.254 

t 

36.00 

C.251 

3.872 

1 

38.00 

-0.069 

2.334 

40.00 

-0.693 

0.319 

] 

42.00 

-1.137 

-0.774 

44.00 

-0.389 

2.218 

i 

46.00 

0.455 

4.772 

48.00 

-0.324 

1.495 

50.00 

-I . 163 

-1.341 

PICKS  OPTION 

ANCLE  OF  ATTACK 
(DECREES)) 

, CY 

CETAMC 

C'’--  2.00FT 

24.00 

0.694 

1.728 

t 

25.00 

0.489 

1.273 

26.00 

0.222 

• .767 

28 . CO 

-0.410 

-0.  36  •> 

30.00 

-1.110 

-1.299 

32.00 

-1.008 

0.441 

34.00 

-0.344 

3.221 

36.00 

0.358 

5.5'  » 

38.00 

-0.099 

3.336 

40.00 

-0.991 

0.456 

42.00 

-1.625 

-1.105 

44.00 

-0.556 

3.169 

46.00 

0 651 

6.819 

48.00 

-0.463 

2.136 

50.00 

-1.662 

-1.916 
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Body  Plus  Tall 
Input : 


VINF 

• 

921.0 

f t/aec 

FSMN 

Mi 

0.8 

RHOINF 

■ 

0.00059 

lb  anc^/ft* 

ANU 

■ 

0,000574 

ft^/sec 

0 

m 

0.312 

ft 

SREF 

m 

0.0768 

ft^ 

Mbs  EL 

m 

0.938 

ft 

BODYt 

m 

3.12 

ft 

XMC 

m 

2.0 

ft 

DELTA 

m 

9.45 

deg 

NS 

m 

50 

lOCONF 

m 

1 

XF 

m 

-2.709 

ft 

YF 

« 

0.0 

ft 

ZF 

m 

0.1561 

ft 

XFMAX 

» 

-2.917 

ft 

YFMAX 

m 

0.0 

ft 

Z^'MAX 

« 

0.312  , 

ft 

SRF.FT 

= 

0.049 

ft^ 

NAOA 

• 

18 

NTYPE 

- 

2 

MIAM 

4 

AO  AD 

rt 

25..  27. 

43.  . 45. 

. 29.,  31.. 

. 47..  49., 

RA 

> 

0.0.  90. 

, 180.,  270 

GAMLIM 

I* 

100  . 

33..  35.. 

51..  53.. 


37. . 39.,  41., 

55.,  57.,  59. 


lOPTl  • 0 

BRN  = 0 

IRHN  . = 0 
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OUTPUT 


ANCLE  Of  ATTACK 
(DECREES) 

trv 

•tfrAm: 

CC-  2.00  FT 

HOt.1.  MilM. 
(FT-U) 

NDKHAL  FORCE 
FIN  MO.  1 (LB) 

2S.00 

0.355 

-0.820 

-0.607 

0.939 

29.00 

0.353 

-1.026 

-0.797 

1.767 

31.00 

0.490 

-0.484 

0.001 

1.120  - 

33.00 

0.767 

0.376 

0.116 

1.208 

35.00 

1.013 

1.171 

0.487 

0.594 

37. C. 

1.190 

1.856 

1.199 

-0.964 

39. Oi 

l.i85 

2.314 

2.215 

-3.388 

43.00 

0.651 

2.669 

-0.155 

-1.884 

45.00 

-0.096 

1.650 

-0.839 

-0.640 

47.00 

-0.804 

1.111 

-2.17,4 

2.498 

51.00 

-1.315 

1.445 

0.374 

1.705 

53.00 

-0.522 

3.404 

2.680 

-0.899 

55.00 

0.176 

4.468 

3.296 

-5.729 

57,00 

-0.238 

3.576 

-0.470 

-1.280 

59.00 

-0.538 

1.324 

-2.087 

2.631 

Th«  results  of 

data  In  Figures 

these  two  sample  outputs,  are  compared  against  experimental 

12a  and  b and  13a  and  b 

6.5  Program  Llml tat  Iona 

Thara  ara  cartain  limits  which  apply  to  the  application  of  the  pro- 
eadurea  and  program  described  in  this  document.  Following  is  a list  of 
the  program  limitations  which  the  user  must  be  familiar  with  before  trying 
to  use  the  program. 

a.  The  program  is  crossflow  Mach  number  limited  to  values  between 
O.IS  and  0>8. 

b.  Due  to  the  large  side  wash  velocities  and  correspondingly  large  fin 
normal  forces  produced  when  a '/o-'tex  core  passes  close  to  a fin«  the  pro- 
gram has  been  set  up  to  discontinue  calculations  at  angles  of  attack  when 
a vortex  core  passes  within  1/4  body  radius  of  a fin  and  proceed  to  the 
next  angle  of  attack.  The  1/4  body  radius  limitation  was  selected  after 
corparing  the  results  of  numerous  runs  against  experimental  data. 

5 potential  vertex  model  is  used  in  the  program.  The  velocity  pro- 
duced by  such  a vortex  at  a radius  P is  given  by  the  following  equation: 

V - r 

2'^r 

It  is  because  of  this  that  the  induced  side  wash  velocities  becomes 
so  large  when  a vortex  core  passes  close  to  a fin.  In  actuality  a vortex 
has  a viscous  core .beyond  which  the  vortex  can  be  modeled  by  a potential 
vertex.  In  the  viscous  core  the  velocity  goes  to  rero  at  the  center.  In 
order  to  retain  as  simple  a model  as  possible,  it  was  decided  to  place  the 
limitation  described  above  on  the  program. 

c.  The  progr^w  should  only  be  applied  to  configurations  which  have  0 « C* 
or  45”.  This  is  because  the  progra.n  rnly  calculates  fin  forces  Induced  by 
the  presence  of  am  asymmetric  vortex  system.  Only  at  these  toll  angles  do 
the  forces  op  the  fins  produced  by  floo  conditions  other  than  vortices  cancel 
themselves  out  due  to  opposite  senses  of  direction. 


d.  Physically  it  la  logical  that  the  side  force  induced  on  the 
finlesa  portion  of  a body  is  contributed  to  by  both  feeding  sheets 
and  shed  vortices.  However,  comparisons  between  predictions  and 
experimental  data  have  Indicated  that  the  best  results  will  be  obtained 
when  the  contributions  from  shed  vortices  are  neglected.  The  program 
logic  predlctr  both  contributions  but  prior  to  the  point  at  which 
total  side  force  is  calculated  the  contributions  from  shed  vortices 
are  set  equal  to  xero.  This  is  accomplished  by  inserting  the  following 


SFPC  (JT)  - 0.0 
FMPGC  (JT)  - 0.0 


Side  force  and  yawing  moment  induced  by  a shed  vortex' 
passing  over  a complete  segment 


SFPSV  (JT)  - 0.0  I 

I Side  force  and  yawing  moment  induced  by  a vortex 
FMGPSV  (JT)  • 0.0  I shedding  within  the  bot  sda  of  a sejment 


within  the  do-loop  ending  with  statement  239. 

These  cards  may  be  removed  at  anytime  in  order  to  study  the  effects  of 
shed  vortices  on  body  side  forces. 

e.  Caution  is  j.dvlsed  when  employing  the  option  to  correct  for 
nose  fineness  and  blunt  ness  ratios.  The  reason  for  this  being  that 
the  scaling  factors  were  derived  from  data  shovrlng  considerable 
scatter.  Also  due  to  the  unsteady  i.ature  of  the  vortex  phenomenon, 
it  is  possible  that  considerably  different  scaling  factors  could  be 
derived  using  other  data  sources. 
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